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Abstract: Using the recently proposed set of discontinuity relations we translate the 
AdS/CFT Y-system to TBA integral equations and quantization conditions for a large 
subset of excited states from the 51(2) sector of the AdS^ x string o"- model. Our derivation 
provides an analytic proof of the fact that the exact Bethe equations reduce to the Beisert- 
Staudacher equations in the asymptotic limit. We also construct the corresponding T- 
system and show that in the language of T-functions the energy formula reduces to a 
single term which depends on a single T-function. 
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1. Introduction 



An important problem in the AdS/CFT correspondence [|| is the calculation of anomalous 
dimensions in the planar A/" = 4 super Yang-Mills theory (SYM) or equivalently the energies 
of the dual AdS^ x string sigma-model. 

The integrability discovered on both sides of the correspondence^ provided us with an 
efficient mathematical apparatus to compute the exact spectrum of the planar AdS / CFT 
models. As the central object of integrability the 2-particle S-matrix of the model ^, |5[ 
plays a prominent role^ and is indispensable for the methods which enable us to compute 
the exact spectrum. 

First, the spectrum of long operators or equivalently states with large R-charge J was 
determined by the asymptotic Bethe Ansatz (ABA) ^, which describes all polynomial 
corrections in 1/J. 

Later, the leading exponentially small corrections^ in J could be taken into account 
by means of the generalized Liischer formulae ||9|, 11]. For the Konishi field their 
small coupling expansion led to beautiful agreement with direct 4-loop field theoretical 



computations |12, 13]. For a certain class of states in the s[(2) sector the 4- and 5-loop 



expansion of these formulae ]14, 15, 16] also satisfy nontrivial consistency checks dictated 
by perturbation theory considerations in the planar = 4 SYM |17, 18]. 



The exact energies, which resum all wrapping corrections in J, can be obtained by the 



application of the Thermodynamic Bethe Ansatz method |1£] to the doubly Wick rotated 
AdS^ X 5^ string sigma-model called the mirror model ]2C, Bl|. 



^For a comprehensive recent collection of review papers, see ]^. 
^For a recent review see |Q, ^ and references therein. 
^These are the so called wrapping corrections. 
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Strictly speaking the TBA method provides only the ground state energy of the model 
in finite volume and its extension to excited states is only a conjecture even if in most of 



the cases it rests on solid grounds. Based on the corresponding string-hypothesis [22| the 
ground state TBA equations of the AdS^ x mirror model were constructed in |23, 24, 25 



and simplified in [26|. In AdS/CFT the ground state TBA equations do not give much 
information about the spectrum since the ground state is protected by supersymmetry, i.e. 
Eo{L) = 0. 

The ground state equations are important nevertheless because they serve as start- 
ing point for the excited state equations. An important further discovery is that the 
Y-functions^ of the mirror TBA equations satisfy the so called Y-system functional equa- 



tions [27|. The Y- functions can be rephrased in the language of T- functions satisfying the 



so-called T-system. The T-system of AdS/CFT lives on a T-hook |27]. The discovery of 
the Y- and T-systems (see Figs. 1, 2) made it possible to determine the asymptotic (large 
J) solutions of the excited state TBA problem. The asymptotic solution is constructed so 
that the ABA and the generalized Liischer formulae are reproduced. 

Based on previously elaborated examples in lattice models ^] and in relativistic 
quantum field theory (QFT) p^, 31, 32, p^], the common experience is that excited state 



TBA equations differ from the ground state ones only in source terms and in quantization 
conditions imposed on objects appearing in the arguments of these terms. These source 
terms can be found by various methods like analytic continuation in some parameter of 



the model ||30|], deforming the integration contour of the ground state equations 
or transforming the Y-system functional relations to integral equations [ 



2|, H, ||. In 



AdS/CFT the analytic continuation [24, Q and contour deformation [34| methods together 
with requiring consistency with the large J asymptotic solution were successfully applied 
to find the excited state TBA equations for certain states of the s[(2) sector of the model. 
In \3t] a general strategy to construct the TBA equations for all states of the model by the 
contour deformation method is outlined. 

Since the TBA equations in AdS/CFT still cannot be derived from first principles it 
is important to test them carefully. In the strong coupling limit it was shown p7| , [3^ that 
the TBA equations reproduce the 1-loop string energies in the semi-classical limit and the 
strong coupling expansion of the energy of the Konishi state fitted from numerical TBA 



computations 39 1 was found to be consistent with direct string theory computations 
1 4^. At weak couphng 5-loop TBA results for the twist-2 states agree |4^, Q with 
those based on the generalized Liischer formula. 

Though the analytic continuation and contour deformation methods provide the TBA 
equations for excited states it is important to derive the TBA equations from the Y-system 
functional relations as well. The main advantage of the Y-system based method is that 
the infinite Y-system can be solved via the T-system and the related T-Q relations, thus 
it opens the way towards the NLIE formulation of the AdS/CFT spectral problem. 



However, until ref. [|4^ appeared, the Y-system based derivation of the TBA equations 
was impossible because the Y-functions are not meromorphic functions like in the relativis- 



*I.e. unknown functions of the TBA equations 
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tic models and have nontrivial discontinuity structure. The main discovery of ref. |44] was 
that this extra difficulty can be overcome for the ground state if the Y-system equations 
are supplemented by appropriate functional relations for the square-root discontinuities of 
the Y- functions. The Y-system equations supplemented by the discontinuity functional 
relations plus some analyticity assumption on the distribution of zeroes and poles of the Y- 
functions are sufficient to transform the Y-system to TBA integral equations. In this spirit 
the ground state TBA equations (including the dressing kernel) were derived 1 44 ] assuming 



that none of the Y- functions have local^ singularities in the entire complex plane. 



It was conjectured in ref. [44| that the form of the discontinuity functional relations is 
state independent. If true, this allows us to derive excited state TBA equations from the 
Y-system as well. In this paper we carry out the Y-system based derivation of the excited 
state TBA equations. We consider states consisting of fundamental particles only, i.e. we 
assume all particle rapidities are real. As a first step we have checked that the conjecture 
is true in the asymptotic limit and found that the asymptotic solutions given in appendix 
C satisfy the discontinuity relations nontrivially. 

Furthermore we show that we do not need to know the "local" analyticity properties 
of the Y-functions in the entire complex plane, it is sufficient to know their behaviour in 



certain regions near the real line. This means that the assumption of ref. |44] concerning 
the "local" singularities of the ground state Y-functions is too restrictive. Indeed, in [^] it 
has been shown numerically that the ground state Y-functions do have local singularities in 
the complex plane. Reconsidering the derivation it turned out that the local singularities 
are arranged in "complexes" (see ref. |Q ) , lie outside the physical strip and do not modify 
the form of the TBA equations. 

Throughout the paper we assume that within certain regions of the complex plane 
the excited state Y-functions are smooth deformations of the asymptotic solution given in 



1 27]. More precisely, we only discuss that part of the parameter space where the exact 
solution has (qualitatively) the same analytic properties as the corresponding asymptotic 
solution. In practice, if the size of the system (J) together with other quantum numbers 
are fixed, this is realized for small enough coupling g. For larger g the TBA equations for 



the given state may undergo phase transitions similarly to what was found in [34, In 
this paper we restrict our attention to the form of the TBA equations valid in the vicinity 
of the asymptotic solution. Then assuming that the solution of the Y-system is found we 
construct the T-functions in an appropriate gauge. These T-functions are also smooth 
deformations of their asymptotic forms. 

The benefit of introducing the T-functions is 2-fold. On the one hand they serve as in- 
ternal variables in terms of which the discontinuity relations and the derivation of the TBA 
equations simplifies drastically. On the other hand T-functions seem to be more fundamen- 
tal objects from the point of view of the AdS/CFT spectral problem. For example we show 
that the complicated TBA energy formula becomes a very simple expression containing the 
single T-function Ti^. This fact indicates that there is an integrable psu{2, 2|4) spin-chain 
in the background such that the Hamiltonian is related to a transfer matrix of the model 



^Here the word "local" means zeroes or poles. 



-4- 



similarly to the cases of lattice models and lattice regularizations |46| of integrable QFTs 
1 47, 48 1 studied previuosly. 



Our final equations agree with previous results of [34]. Our derivation is based on 
functional relations and analyticity assumptions which we know are satisfied exactly in the 
asymptotic limit. Hence as a by-product our results prove analytically that the asymptotic 
solution satisfies the large J limit of the TEA equations and also that the large J limit of 
the exact Bethe equations coincide with the ABA equations. This fact, although it played 
an important role in the extraction of the 5-loop Liischer terms from TBA |41, 43|, has 
not been proven analytically so far. 

The organization of the paper is as follows. In the next section we introduce our starting 
relations and assumptions. Section 3 contains important lemmas and the transformation of 
the Y-system to TBA integral equations leaving the necessary discontinuities temporarily 
unspecified. In section 4 we construct the T-system in a partially fixed special gauge which 
turns out to be very useful to simplify the derivation of the TBA equations. In section 5 
and 6 we compute the so far undetermined discontinuities y^"^/y^") and A from dispersion 
relations and we fix completely our gauge choice for the T-functions. In section 7 we derive 
the simplified version of the TBA equations, while section 8 contains their canonical and 
hybrid form. In section 9 we discuss the quantization conditions and the exact Bethe 
equations. Finally, in section 10 we discuss a simplified energy formula. The paper is 
closed by our conclusions. Appendix A contains the definitions of frequently used objects 
and kernels, appendix B the discussion of the branch cut discontinuities of the dressing 
part of the discontinuity function A. In appendix C the asymptotic solution of the Y- and 
T-system equations are presented, in appendix D the precise definition of the discontinuity 
relations is given and finally in appendix E we discuss the meaning of the exact Bethe 
quantization conditions. 

2. Starting relations and assumptions 

Our starting point in this paper is the hypothesis that the same Y-system describes all 



excited states in the planar AdS/CFT spectral problem |27]. The Y-system of AdS/CFT 
takes the standard form 

(i + i/y,_i,,)(i + i/y,+i,,)' ^ • ' 



but it lives on an (a, s) lattice represented in Fig. 1. This is equivalent to imposing 
the boundary conditions Yo.s — ^ oo> ^2,|s|>2 ~^ oo and Ya>2,±2 while the product 
^,±2 ^,±3 should be kept finite in order to l2,±2 be finite. This last requirement shows 
that the Y-system defined in the domain of Fig. 1 cannot form a closed set of equations 
because the new functions l3,±2 y2,±3 enter the problem and we need additional equations 
(independent of the Y-system) to complete the system. 

The Y-system is equivalent to a T-system defined on a T-hook of Fig. 2: 

'^a.s '^a,s — Ta+l,s Ta-l,s + 7a,s+l 7a,s-l, (2.2) 
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where the relation to the Ya.s functions is given by 

rp rp T-'+ rp— rp-\- rp — 

-ta,s+l -ta,s— 1 1 I V a,s a,s 1 _l 1 /V a,s a,s q\ 

a,s — 7^ Tf, ) + ^ a,s — Tf, Tf, ) -'- + -'-/ ^a,s — Tf, Tf, • I'^-'^J 

^ a+l,s J- a—l,s ^ a+l,s a—l,s ^ a,s~l a,s+l 

The T-equations can be extended to the infinite (a, s) lattice by imposing the boundary 
conditions that all Ta^s = outside the T-hook. The T-system is a more fundamental set of 
equations than the Y-system because it forms a closed set of functional equations^ which 
determines all the Y-functions together with the supplementary combination 5^,±2^2,±3- 
The T-equations and the Y-functions given by (^]^) are invariant with respect to the gauge 
transformations'' 

Ta,s ^ 5^^' gf^ gf^ gi-""^ (2.4) 

with (?i,2,3,4 being arbitrary functions. In this paper we will choose a gauge where Tq^s = 1. 
This fixes gigs = g2gi = 1. 

It is known that in AdS/CFT the Y-functions have square root branch cuts and live 



on an infinite genus Riemann surface |25, 26, 44, 49 1. The different sheets of the Riemann 



surface are connected through square-root branch cuts starting from the branch points with 
real parts ±2 and run to infinity along the lines with integer (in 1/g units) imaginary parts. 
To avoid the complications connected to using different sheets we will use a convention 
where all our functions are assumed to be defined on the first Riemann sheet, defined as 
the entire complex plane excluding the above cuts. In our convention if we analytically 
continue a function through one of the cuts, this becomes a new function, which can also be 
continued to the entire first Riemann sheet. We will often use the operation f{u) — )• f*{u) 
corresponding to analytic continuation of the function f{u) through the real cut \u\ > 2. 
The construction of f*{u) consists of the two steps of analytic continuation of f{u) through 
the real cut followed by analytic extension of the new function to the entire first Riemann 
sheet. Our functions are assumed to be meromorphic in the first Riemann sheet: they 
have discontinuities along some (but not necessarily all) of the cuts of the first Riemann 
sheet and cannot have other discontinuities but may have local singularities (zeroes and/or 
poles) in the sheet (and in some cases even on the cuts). 

Restricting the Y-functions to the first Riemann sheet the Y-system equations supple- 
mented by some analyticity information on the local singularities are not enough to derive 
the TBA integral equations. Some additional information on the discontinuities is needed 
as well. According to the proposal of ref. |44] this missing piece of information is a set of 



functional equations relating the discontinuities in a state independent way. These discon- 
tinuity equations translated to dispersion relations determine the extra 13^2 ^2,3 functions 
as well. 

We will use the conventions of ref. [ ]34| throughout the paper. The Y-functions in 
these conventions are related to the y^.s variables as follows: 

yi"^ = Yl"^ = -Y2,a2, a = ± Yq = Yq,o Q = 1,2, ... (2.5) 



^Here the word "closed" means that there are as many functional equations as T-functions. 
''Our notations and conventions are explained in appendix A. 
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Y 



(a) 



l/Yr 



Y 



(a) 



Y 



a 



± 



m 



1,2,... 



(2.6) 



'm\vw ~ "-I ^■m+L,al ^ m\w ~ l,a(m+l) 

From the ground state equations we obtain the following structure for the locations of the 
branch cuts: 



Y^\u) : u + i2m/g m £ 



Ym(u) 



u±i{m + 2j)/g, 



j = 1,2, ...,oo, 



where u G (— (X),— 2)U(2,oo). One of our main assumptions is that this structure remains 
valid for the excited states as well. For later convenience we rewrite the discontinuity 
relations proposed in for the square-root branch cuts in the conventions we are using 
in this paper. We introduce the symbol [f]z with Z £ Z to denote the discontinuity of the 
function f{u): 



[fin)]z = 
and define 



lim {f{u + i Z/g + ie) - f{u + i Z/g - ie)) , 

E^0 + 



[inri 



lJ+1- 



u £ (-00,-2) U (2, do) (2.7) 



(2. 



The discontinuity relations relate the "jumps" of this function and those of some other 
Y- functions. They take the form 



[A] 



±2N 



In 1 



Y- 



(a) 



N 



±2N 



m=l 



In 1 + 



1 



Y 



(a) 




\ 



I 



In 



Y 



Y_ 



N 



±(2Af-Q) 



(2.9) 
(2.10) 



±2N 



with iV = 1,2, ... and 



lluui 



In 



l-Y 



(«)' 



±1 



l-Y 



(a) 



InK 



(a) 



l\w 



In 



±1 



'i-i/yi")' 



(2.11) 



In the sequel [/(ii)]z means the analytic extension of the discontinuity (|2.7| ) to generic 
values of u. Using the notation defined earlier in this section we can write [f{u)]z = 
/[+^l(n) — {f^^^^)*{u) where (/'"*" ^')*(n) is the function obtained by analytic continuation 
of (/'^^')('u) through the cut lying along the real line (crossing it from below) and then 
extended to the first Riemann sheet. 

Some important remarks on the interpretation of the discontinuity relations (^-2.11) 
are in order. As they stand, (2.8- p.li ) are valid if the Y-function combinations appearing 
in them have no logarithmic discontinuities crossing the lines of the square-root disconti- 
nuities. In order to get rid of the difficulties caused by the logarithmic discontinuities we 
use the derivative of the discontinuity relations (p. ^ 2. 11 ) as our starting point and apply 
dispersion relations for the derivatives of A and lnY_^^ /Yj"\ Then a second subtlety 
appears if there are local singularities of the Y-function combinations lying exactly on the 
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lines of square-root discontinuities. Such local singularities do not modify the disconti- 
nuity relations, but they contribute to the corresponding dispersion relation through the 
residue theorem. To cure this problem a finer interpretation of the discontinuity relations 
is necessary in which the contribution of such local singularities are taken into account as 
well. 

To find the correct interpretation we can invoke the asymptotic solution. It can be 
shown that only ( p.9[ ) at = 1 and at the positions of real poles of 1 — l/y_^"^ must 



In 



1 + l/Yl 



(a) 
\vw 



±1 



on the right hand 



be refined. The term which causes the trouble is 

side of ( |2.9D , because its derivative has poles sitting right on the discontinuities^ with 
Imu = ±l/g(. This refined interpretation of (|2.8 - p.H ) is necessary only when they are 
translated to dispersion relations. In order for the dispersion relation applied for the 
derivative of ( |2.9| ) give the correct formula for A the following replacement must be done 
on the right hand side of (|2.9|) at = 1: 



In 1 + 



1 



y/,"^ (nl 



±1 



In 



1 + 



1 



1 



y/°) (n) 



(2.12) 



J ±1 



where {u) is the polynomial having zeroes at the positions of the real zeroes of y^"^ 
with absolute values larger than 2. In other words the poles corresponding to the zeroes 
of P2^\u) must be ignored. For a proof of this formula for the most general state of the 
model see appendix D. 

We note that as a consequence of the definition ( |2.8D , and the fact that going around 
twice a square root branch point gives back the original function, we have the relation 



A(n) = -A^u). 



(2.13) 



The Y-system (2.1) and the discontinuity relations (^- 2.11] ) are not sufficient to derive 
the TBA equations. Some more information is needed about the discontinuities lying along 
the real axis. Based on the properties of the solution for the ground state TBA equations 
and the similar properties of the asymptotic solution for excited states, we require that the 
"fermionic" Y- functions are analytic continuations of each other: 



(2.14) 



We will also assume that all Y- functions are real analytic. This assumption is based 
on the observations that the Y-functions are real analytic functions for the ground state 
solution of the TBA equations and also in the asymptotic limit for excited state solutions. 
Since we consider the AdS/CFT Y-functions as smooth deformations of their asymptotic 
counterparts we restrict ourselves to deformations that preserve the property of real ana- 
lyticity. We note that in the absence of this assumption several TBA integral equations 
could be set up (such that their asymptotic limit coincides). 

To summarize we assume that 



*We note that other Y-combinations for other values of N can also have local singularities along the 
cuts under consideration, but they mutually cancel each other's contribution. 
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• Y- functions satisfy the Y-system equations (|2.l| ). 

• The discontinuities are related by the modified relations (2.8- 2. 11 ), ( 2.12| ). 

• Yi"'\u) and Yi°'\u) are related by 

• The Y- functions are real analytic functions. 

• The Y- functions are smooth deformations of their asymptotic limit. 

If one assumes that both A and y^'^^/y^'^) are analytic and bounded near ±2 then their 
discontinuities must be zero at the branching points ±2. From this requirement and from 
( |2l3D , ( ^l4D it follows that the combinations A{u)/V^-u'^ and ln(yi"V^|"V\/4 - 



have no square root branch cuts along the real axis. These properties have been used in 



the derivation given in ref. [44| and we think it is important to emphasize them since they 
play very important role also in our considerations. 

The form of the discontinuity relations (2.8- p.ll ) and ( p. 14 ) have been conjectured 
to be independent of the particular excited state of the AdS^ x sigma-model under 
consideration. Using the formulae and results of appendices A, B, and C this conjecture 
can be proven to be valid for the asymptotic solutions. Here we will verify^ the relations 
( 2. 10 ), ( 2. 11 ) and ( 2.14|) in the asymptotic limit The (asymptotic) justification of ( ^^ ) 
will be presented in section 6. 

With the help of the T-function representation of y!^^^ /Y^^ and the formulae of ap- 
pendix C it can be shown that asymptotically: 

infyi°V>^f ^ 



In "^"^ 



R 



R 

m -Dp 



(2.15) 



from which 
immediately. 

Next we write 



l^(y(0)/y(0) 



±2N 



0, the asymptotic (Yq — )• 0) limit of (|2.10| ), follows 



1 



1 



Y_ 



(0) 



7^(0) 
^2,1 



^{0)+^(0)- 
-^2,2 -^2,2 
7.(0)7.(0) 
2,1 2,3 



1 + 



(0) 



2,2 



Y 



(0) 



(2.16) 



where we have used the identity 



7.(0)+ 7.(0)- _ T^(0)+ ^(0-7.(0) 

-^2,2 -^2,2 — ^ ^2,3 



(2.17) 



and this proves ( |2.14 ) asymptotically. 



The asymptotic verification of (2.11) goes as follows. Using the T-representation: 



InY 



(0) 



l\vw 



±1 



Iny 



(0) 



2,1 



±1 



In 



7.(0) 
-^1,1 

7.(0) 

-^2,2 J ±1 



(2.18) 



^During the verification the upper case index (a) is ignored since in the limit under consideration 
the two wings of the Y-system become independent. Furthermore we forget about possible logarithmic 



discontinuities since (2.8-2.11) account for the contribution of the square-root discontinuities only. 
^''Asymptotic limit: J — 00 or g 



0. 
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With the help of (|a6|)- (|C.14D it can be shown that 



In 



r(o) 

4,2 



In 



±1 



^(0)+ ^(0)- ^(0) 
-^1,1 -^1,1 -^3,2 

Jc-(0)+ gio)- ^ 



In 



1 + 



In- 



1 -y 



(0) 



1 -y 



(0)- 



(2.19) 



Formulae ( p.l8 ) and ( 2.19D together verify the first relation of (2.11) while the second 
relation follows from ( 2.19| ) and the following two relations: 



iny 



(0) 



1 to 



±1 



Iny 



(0) 



1,2 



±1 



In 



7^(0) 
-^2,2 



+ 



InT, 



(0) 
1,3 



±1 



±1 



(2.20) 



InT, 



(0) 
1,3 



±1 



-'I,! -'2,2 



In- 



y 



(0) 



y 



(2.21) 



In the derivation of the excited state TBA equations we can avoid the technical prob- 
lems coming from dealing with branch cuts of the log functions if we first derive equations 
for the derivatives of the logy-functions. To this end we use the logarithmic derivative of 
the Y-system (2.1), and the derivative of the discontinuity relations (^l8- 2.1l] ), ( |2.13| ) and 
(2.14). Qualitative information about the local singularities of the Y- functions can be read 
off from their asymptotic form. Finally we integrate the equations for the derivatives. 

In order to be able to extract the necessary analyticity information on the local singu- 
larities of the Y-functions we will make use of the assumption that the exact Y-functions 
are smooth deformations of the asymptotic ones. This however cannot be satisfied in the 
whole complex plane, but only in certain strips (similarity regions). We will show that 
analyticity information regarding the behavior of our functions in these strips are enough 
to derive the TBA equations and to determine the exact Y-functions. 

The Y-functions are smooth deformations of their asymptotic counterparts as long as 
the y^°^ funct ions are small. From (|C.26| ) it follows that this condition is satisfied in the 

region = {u G C : llmuj < Q/g, \u — i^f^^^ > eg}, where eg is a small but 
not infinitesimal positive parameter for Q > 2 and zero for Q = 1. Then the Y-system 
equations (2.1) imply for the other Y-functions the following similarity regions: 



Y 



(a) 



m\vw 



(n) 



yi"Vn) 



y(^l(°^(n) 



y 



' m\w 
(a)(0) 



y 



(0) 



l,Q:(m+l)('"') 



liGl^eg, (5 = 1,2,..., 

Ilmul < mjg^ 
Imn| < mjg^ 
|Imu| < 2/<^, 



m 



1,2,..., 



a 



1,2, 



yf)(0)(n) = -y2^;^yn), |Imn|<2/5, 



.(0) 



a 



The discontinuity relations (|2.8| - p.ll ) simplify considerably in the language of T- 
functions. The complicated products of Y-functions in the argument of the log function 
become a product of a few T-functions only. In section 4 we will show that there is a 
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Figure 1: The AdSs/CFT^^ Y-system. Full circles on the a axis correspond to massive nodes with s = 0. 
The s axis corresponds to nodes with a = 1. 




Figure 2: The AdSs/CFT^ T-system. The s axis corresponds to nodes with a = Q. 



particular choice for the gauge where the exact T-functions are smooth deformations of the 
asymptotic ones and due to their simple square-root branch cut structure the right hand 
sides of (2.8-2.11) simplify drastically. 



3. TBA equations with cuts 

In this section we transform the Y-system equations ( p.l| ) into TBA integral equations. 
This transformation is not as complete here as for the case of integrable relativistic models 
because of the presence of cuts (discontinuities) in the analytic extension of some of the 
Y- functions. The Y-system equations are of the universal form 

y+y-=n, (3.1) 

but the details of the corresponding integral equation depend on the analytic properties of 
the "unknown" function y. 
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3.1 TBA lemma la 



Assume that the set of zeroes [poles] of y{u) mside the physical strip {—1/g < Imu < 1/g) 
is {Cj} [{%}]• We also assume that there may be discontinuities along the cuts with 
imaginary part |: 

y(u + - + ie) . , 

? = e^("+"\ \u\ > 2, V{±2) = 0. (3.2) 

y{u + ^- ie) 

Further we assume the large u asymptotic behaviour 

y{u) you^'°, u — ^ oo (3.3) 

below the cut and 

y{u) ^ y+ n*^+ , u — > oo + - + ie (3.4) 

9 

just above the cut. (If there is no discontinuity then of course yo = ?/+, Mq = M+.) The 



Y-system equation (3.1) implies that Tl{u) behaves asymptotically as 

TZ{u) « y+?/on^^"+^+ , u — ^ oo + ie. (3.5) 

The function 

t{u) = , 7, (3.6) 



where 



t(n)=tanh^ (3.7) 



has the same zeroes [poles] (in the physical strip) as y{u) and satisfies r^r = 1. Using 
this function and the universal TBA kernel function 



the solution of (3J) which has the analytic properties detailed above is 

y = sgn(yo)Texp{ln7l'*s - * s}. (3.9) 

Here TZ^{u) = TZ{u + ie) and V''{u) = V{u + ie) and this notation indicates that the 
integration contour in (|3.9| ) goes just above the real line. Our definition of the log (In) 
function is the standard one: we always assume that the cut of this function is along the 
negative real axis. 

The result (|3.9| ) can be proven directly by substituting this expression into the func- 
tional equation ( |3.lD . 
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3.2 TBA lemma lb 

In one of the TBA equations the Y-function has no cuts near the physical strip but singular 
points right at its boundaries. In this case we assume again that the set of zeroes [poles] 
of y{u) inside the physical strip is {^j} [{%}] and further assume that the set of zeroes 
[poles] on the boundary of the physical strip is {za ± ^} [{wp ± ^}], where Za [w/?] are real. 
The large u asymptotic behaviour is 

y{u)^you^'°, u — ^ oo. (3.10) 



The Y-system equation (|3.1|) implies that TZ{u) has double zeroes [poles] at {za} [{wp}] 



and behaves asymptotically as 

n{u) ^ ylu"^^^" , u — ^oo + ie. (3.11) 
In this case the solution of (|3.1| ) with the right analytic properties is 

y = sgn{yo)T exp{lnTZ -k s} . (3.12) 



The result (p^ ) can be proven directly because the logarithmic singularities (at {za} 



[{wis}] or at oo) are integrable and do not invalidate the result. 
3.3 TBA lemma 2 

Let us now assume that the "unknown" function y{u) has zeroes [poles] inside the physical 
strip as before, and in addition has (multiplicative) discontinuities along the real cuts: 

y{u + ie) ^^Jiu+^e)^ |^| > 3, J(±2) = 0. (3.13) 

y{u - le) 

We also assume y{u) has constant asymptotics for large u: it approaches the constants 
y+[y_] when u — > oo just above [below] the real line. This implies 

e^("+-) « Jo, = ^, 7l(n)«y+y_, u^oo. (3.14) 

y- 



In this case the solution of ( |3.1| ) is 

y = sgn(y+)r expjlnT?, * s + * si}, (3.15) 

where 

= 4smh(S-^e) - 

3.4 TBA integral equations 

Using the above two lemmas we now write down the set of TBA integral equations corre- 
sponding to the Y-system ( |2.lD . First we spell out these equations using the notations of 
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mlvw m\vw 



l|t;ui l\vw 



l|ui l|ui 



(1 + Y^+i) 
1 _ yj") (1 + Y2) ' 



m > 2, 



m > 2, 



y(") 
1 _ ^ 



y{+) y(-) 



-1 



(i + yQ+i)(i + yo 



1 + r^^"^ 

(i + yi)(i + y/i:,Y 



(3.17) 

(3.18) 
(3.19) 

(3.20) 

Q > 2, (3.21) 
(3.22) 
(3.23) 



Using the notation X for a general index, which can take the values n, , 

or l^'*, and denoting the set of zeroes of Yx in the physical strip by and the set of 

its poles by {r/x.fc}, finally the sign of Yx{u) in the limit u — > cxd + ie by (sgn)x, we define 



tx{u) = (sgn)x 



Y\,t{u-U,j) \ 
Wk - Vx,k) J ■ 



(3.24) 



From the discontinuity relations we see that we can use Lemma la with V = for the 
cases ( pTl) , (p^ ) and if Q > 3 also for (^^). Lemma lb is used for (^^) if Q = 2. 
Further we have to use Lemma la 



for the ^1^''"^ equation (3.18) with V = In 



for the Y^^ equation (|3.20[) with y = In I - 




for the Yi equation lK2^ with V = A 



and finally Lemma 2 for the equation ( |3.23| ) with J = 
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We find 



Y 



(a) _ Aa) 



m\vw m\vw 



.(a) _ Ja) 



t I exp < In 



^ m+l|'!;«;' ^ m—l\vw' 



(a) 



(1 + Wl) 



* S 



m>2, 



* s + In 



1 - Y_ 



(a) 



1 - Y_ 



(a) 



■k S 



y{a) _ (a) J In [l I 



* s + In 



"1 - 


1 -] 

y(") 




1 - 


1 


v(q) 






YQ+,YQ_,{l + Y^t\U{l + Y^-_\ 



Q — l\vw' 



^2 (i-Fi+^)(i-Fro 



y(+) y(-) 



★ S — A -Jk- S 



(i + ri)(i + y/i:J 




(3.25) 

(3.26) 
(3.27) 
(3.28) 

Q>2, (3.29) 
(3.30) 
(3.31) 



There is no equation for 1^"'', but since it is the analytic continuation of Y!f\ the 
very definition of the discontinuity J^°) can be written as 

^y(a)e^p|_j(a)| (3 32) 

The above set of TBA integral equations is incomplete yet, since the discontinuities A and 
are undetermined. They will be obtained using dispersion relations in Section 6 and 
Section 5, respectively. 

4. Construction of the T-system 

In this section we start to construct a T-system from the Y-system. This is not unique 
because there is a gauge freedom in this transformation. Since our main assumption is 
that the exact Y-functions are (at least qualitatively) close to their asymptotic values 
given by the Bethe Ansatz solution and in our subsequent considerations the T-system 
plays an important role, we will choose a gauge in which also the T-system is close to the 
asymptotic solution. This will be achieved here only partially and the final, complete gauge 
fixing will be given in Section 6. 

4.1 Chain lemma 

Let us assume that we want to find the solution of the infinite system 



a = 1,2,..., 0-0 = 1, 



(4.1) 
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where the unknowns (Ta{u) are assumed to have no zeroes/poles near the physical strip (in 
the physical strip and a little beyond) . (u) on the right hand side are given functions not 
having zeroes/poles near the real axis and taking complex values (excluding the negative 
real axis). Furthermore, (only) may have discontinuities on the real axis: 

6(n + ie) _ j(^u+ie) 



^i(u - ie) 



e'^™, /(±2)=0. (4.2) 



In such a case we assume that (only) o"i has discontinuities along the cuts with imaginary 
part 

cji(n + - + ze) / ,.N (Ti(u— - + ie) , ,■^ 

cri{u + l^-ie) ' cri(n - ^ - ze) 

where 

p + m = I, p(±2) = m(±2) = 0. (4.4) 



The solution of ( [4.1[) is given by 

a, = exp|f;ine^*£^-p^*i^„| (4.5) 

where 

A-l 

e^=Y,Ka+l-A+2j. (4.6) 
j=0 

(Here we use the convention K-a = —Ka, Kq = 0.) Note that 

£^=e'\, £l = Ka. (4.7) 

Because of the discontinuities, we have to integrate slightly above the real axis as indicated 

by 

eAiu) =U{u + ie), p'{u) =p{u + ie), (4.8) 

but this is only necessary for ^i. 

4.2 Constructing q 

To construct Tafi we can use the chain lemma with 

(Ta=Ta,0, ^a = l + Ya, = 1,2,... (4.9) 

(We use the Tq^s = 1 gauge throughout this paper.) 

From the asymptotic solution (and assuming similar behavior for the exact solution) 
we see that for a > 2 1 + has no zeroes/poles in the strip with imaginary part in the 
interval (^^i ^^) (the strip (1 — a, a — 1) for short) and no cuts in the strip {—a, a). 

The first singularities are poles at and the first cuts are at it— . 1 + Yi has no 

zeroes/poles/cuts in (—1,1). 

^1 has no real cuts in this case and we can put p = 0. In the first step we construct 
Ta,o by using the chain lemma formula ( [4.51) . This gives the solution along the real axis 
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but it is easily extended to the whole physical strip (and a little beyond). The next step is 
to use the defining relation 

l + Ya, a = 1,2,..., (4.10) 



ta+l.O-ta-l.O 



which (for a > 2) extends the solution so that Tafi is free of any zeroes/poles in the strip 
(—a, a) and meromorphic in (—1 — a, a + 1). (The first poles are at i^^"^ and the first cuts 



are at ±ii^^ili.) Ti^q is more regular: it has no zeroes/poles/cuts in the strip (—2,2) with 
first cuts at ±y. 

4.3 Constructing Tq i 

Here we restrict our attention to the s[(2) (sub-)sector (see appendix C) only. In this 
special case the two sides of the Y-system are identical: 

^^1 = ^-1--' ^Hi = ^H-' Y±'^=Y±^ m = l,2,..., a = ±. (4.11) 
We can use here the chain lemma using the identifications 

Ta,i = -aaTa{-l)\ a = 1,2,..., Tq,! = 1 (4.12) 

and 

ia = Ta+lTa-l{l + y^^, a = 2, 3, . . . , = -T2 (^1 - (4.13) 

In this case there are real cuts for ,^1 and to specify the solution completely we also 
need p{u). This function will be fixed later (in section 6). Again, we first use the formula 
(|4.5| ) to determine 1 near the physical strip with Ti 1 having discontinuities along the 
=b^ cuts. Then, using the defining relation 

l + y,,i, a = 1,2,..., (4.14) 



fa+l.l^a-l,! 



we can extend the solution. The Ta^i constructed this way will be meromorphic in the strip 
(— o, o) and the first cuts occur at it—. 

V ) / g 

4.4 The complete T-system 

Having constructed the T-system elements Tafi and T^^i, the rest of the T-system can 
simply be calculated from the relation between the T-system and Y-system elements. For 
example, we have 

Ta,2 = T? Tf, , a = 2,3,... (4.15) 

^ a—l\vw ^a,0 

From this representation we can see that Ta^2 is meromorphic in the strip (1 — a, a — 1). In 

the a = 1 case we have ^ 

To 1 1 

T^,2 = -T^Tr- (4.16) 



-17- 



This function has discontinuities along the real cuts inherited from YL. The other factor 
is a meromorphic function in the (—2,2) strip. 

It is also possible to calculate Ti^g, T2^s for s = 3,4,... and Ta.s for s < 0. These 
functions will not be used in our considerations. We just note that Ta^s / Ta-s in general, 
in spite of the fact that the two sides of the Y-system are identical in the s[(2) (sub-)sector 
we are considering here. To illustrate this, we calculate 

Ta-i = a = 1,2,... (4.17) 

from the relation between the T-system and Y-system elements. Although it is not at 
all obvious from the above formula, we know that by construction the set {Ta-i} must 



also satisfy ( |4.14 ) with the same right hand side (assuming that we stay in the 5l(2) 
(sub-)sector). This means that the T-functions {Ta-i} are gauge transforms of the func- 
tions {Ta^i} in the sense discussed below. Indeed, using the results given in appendix 



C, we can verify the above structure by explicitly calculating {Ta-i\ from (|4.17D in the 
asymptotic limit. We also find that the T-functions {Ta-i} are exponentially small asymp- 
totically. 

4.5 Gauge transformations 

The relation between the Y-system and T-system is not unique, there is a gauge freedom 
Ta^s — ^ jfl.s- We restrict this gauge freedom by demanding 



1, Ta,o = Ta,o, (4.18) 



i.e. we work in the Tq s = 1 gauge and also fix Ta^ as constructed explicitly above using 
the chain lemma. The remaining gauge freedom is of the form 

^'^''^ = f[a+s] fl-a-s] ^^-^^^ 

The gauge transformation is generated by a single function /3(u): 

ri,i = /3ri,i, p = —^—, (4.20) 

T2,i =/3+/3-r2,i, r2,2 = /3++/32/3-~r2,2, ... (4.21) 

We see that the Ta^i solution constructed in this section with some p{u) is of the form 
as if it were a gauge transform of the p{u) = solution with the discontinuous gauge 
transformation 

^ = ^1. ,(„,=„pUi^:i*lj. (4.22) 

q \^J ZTTl u — V ) 

Since our Y-functions and T-functions have discontinuities, it is natural to allow also gauge 
transformations that have discontinuities along the same cut lines. 
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4.6 Large |n| asymptotics 

It is easy to show that if the function f{u) has large \u\ asymptotics 

/(n) /-t + Kin u — > ±00, (4.23) 

then this is reproduced by the / * Ka convolution: 

{f*Ka){u)^f± + Kln\u\, u^±oo, (4.24) 

and similarly for the modified convolution / i Ka- From the known asymptotics of 

and Y- we can determine the large \u\ behavior of the functions on the right hand side of 

( |4.1| ) for the T^,! problem: 



a? B 
^a{u) ^ ^ — -, a = 2,3,..., \u\ — ;> 00 (4.25) 

— 1 u 



and the discontinuity behaves as 

Bu 

p{u) ^ In— . (4.26) 

Thus we have ^ ^ 

In^i^Ka^ln-, pi Ka ^In-^ (4.27) 

u \A\ 

and further 

. . \A\ i ^ , , fm + l)2 1 1\A\a , 
Oa{n) « ^ exp <^ V mm m + 1, a In ^ \ = 4.28 

\m=l ) 

which reproduces the expected large |n| behavior of T^^i. 
5. Dispersion relation for y_ 

In this section we determine the discontinuity J^°^' = In using the discontinuity relation 



( 2.10 ), which can be rewritten as 

N 

±2N 



Here we introduced the notation Lq = ln(l + yQ). Using the notation tg = IuTq^q (^o = 0) 
and the relation 

-^Q = *Q + *Q ~ ~ <3 = 1;2, (5.2) 

( |5.1| ) can be drastically simplified. In the language of the tq variables most of the terms 
cancel and we are left with 



±2Af 



[tN+l]±N ~ [*l]±2Af ~ [*Af]±(Ar_i) - - [tl]±2N ' (^-3) 
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Im 



[4] 



[2] 



Re 



[-4] 



Figure 3: The contour To- It goes around all positive and negative even cuts. 



The last equality follows from using the analytic properties of the Tq^ elements con- 
structed in the previous section by the chain lemma, namely that Tq^ are analytic in the 



Y 



strip {—Q,Q)- The last relation is equivalent to saying that the combination ^-j^Ti^q has 

discontinuities along the real axis only (in the gauge we are using) or equivalently that 
H^'^^ has discontinuities only on the real axis, where 



?2,3 



Hi-) 



T2, 



T3, 



(5.4) 



5.1 Dispersion relation 

Note that 



u + ie) 



\u\ > 2, 



(5.5) 



which follows from the fact that (i±)* = Kp, i. e. Y± are analytic continuations of each 
other. Moreover, since we know that exp{ J^°)} has no real zeroes/poles (this is true in the 
asymptotic limit and according to our basic assumption it remains true also exactly) the 
combination 

J(")'(n) 



V4" 



(5.6) 



is meromorphic near the real line (in the strip |Imti| < ^, where 7 is small, but not 



infinitesimal) and has no other poles than those at n = ±2. Note that the function (| 
goes like l/u^ for large which is necessary for some of our integrals to converge. The 
above properties allow us to define (for < Imn < ^) 



r(")( 



u]= (t> dvji'^)'{v)K{v,u), 

To 



(5.7) 



where the contour Fq goes around all the even cuts ib2A^ (see Fig. |3|) 
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For later use we define also the contour integral z oo]^^^ where the contour consists 
of a horizontal line f + ^ and goes around all the Z + 2m cuts (m = 0, 1, . . . ). Here ■j < Z, 
Z a positive integer. See Fig. ^. The contour integral _z is defined similarly. 

Here the contour goes around all the —{Z + 2m) cuts (m = 0,1,...) and comes back (from 
right to left) along the horizontal line w — y. Here —Z < —7, Z a positive integer. 



Evaluating the integral (|5.7| ), it can be written as a sum of three terms, the first being 
the contribution of the narrow strip — ^ < Imf < ^: 

r(")(n) = (Pj("))('u)+ / dvJ^'^^'{v)K{v,u)+(f dvJ^''^'{v)K{v,u). (5.8) 

J[7,2,oo] J[-oo,-2,-7] 

Here we introduced the derivative operator D defined as 

(Pg)(n) = q'{u) + 77=^= {(u - 2)grf(-2) - {u + 2)qd{2)} , (5.9) 



where qd{u) = V4 — v?q'{u). The above three terms correspond to the three poles [v = u, 
V = ±2) in the narrow strip, assuming that q{u) has square root cuts for \u\ > 2, in which 
case qd{u) is analytic around the branch points u = ±2. 

The usefulness of this definition can be seen from the fact that the T> derivative can 
be "integrated" in the following sense. If the function g{u) can be represented as 

/oo 
dvfiv)K^-^{v,u), (5.10) 
-00 

then 

/oo 
dvf'iv)K^''\v,u). (5.11) 
-00 

If 

h(v) = 

2m 



h{u) = ^.a{u,i), (5.12) 



where ^ is constant, then 

and finally if 

then 



{Vh){u) = -K{i,u) (5.13) 
^(u) = 261n(-x(n)), (5.14) 



{VA){u) = Q, (5.15) 

i. e. ( ^.14| ) is the zero mode of the derivative V. Note that e'^*-^'^^ = 1 if 6 is integer and 
note also that 

^(n) ~ — 261n |u| + const., u — ie it 00. (5.16) 

This last property makes it possible to adjust the large \u\ behavior of the inverse of the 
operator P, by adding a multiple of the zero mode, if necessary. 
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There exists an alternative way of calculating { f).7\) using the (derivative of the) dis- 
continuity relations ( |5.lD . The contribution of the integrals along the cuts in the upper 
half plane is 



dvK^^\v,u)L'Q{v) - i dvK^^\v,u)L'Q{v) 

> ^[0,Q,oo] 



(5.17) 



Here the pole term {(f part) simplifies drastically in terms of tq (using the same analyticity 
information as was already used above) and becomes 



[0,2,00] 



dvK{v,u)t[{v). 



(5.18) 



After similar considerations concerning the contribution of the cuts in the lower half plane, 
we arrive at 



00 n QQ 

r(°)(n) = ^ / dv L'q{v) (k^Q\v,u) - K^-'^\v,u)\ 



(5.19) 



dv t[ {v)K{v, u) - <f) dv t[ {v)K(v, u) . 

'[7,2,00] ^[-00,-2,-7] 

Comparing ( |5.8|) and ( |5.19| ) we can express as 

(Pj("))(n) = / dvL'Q{v)[K^'^\v,u)-K^-^\v,u)\ 

+ * dv In' H^'^\v)K{v,u) + (f dv In' H^'^\v)K{v,u) 

/[7,2,oo] /[-co, -2, -7] 



(5.20) 



We already know that H^"'^ has no discontinuities in the upper and lower half planes, but to 
be able to "integrate" ( |5.2C1| ) we need to know the position of its zeroes/poles. We assume 
that its set of zeroes is {xi"^} and the set of poles is {y^"^} (and also assume that these 
singular points are not real and are not on any of the even cuts). Then, using ( 5.10| - 5?I^ ), 
the "integral" of ( 5.20D can be written as 



J(")(n) = ^a(n,yi"))-5]a(..,xi")) 

b a 

roo p 

+ ^ dvLq^v) K^^\v,u) - K^'^\v,u) 



(5.21) 



Since for large |n| J^"^ — > const., there was not necessary to add the zero mode (|5.14|) to 
this solution. 

According to our basic assumption, H^""^ is close to its asymptotic counterpart: 



(5.22) 
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Im 



IZ+4J 



[Z+2] 



[Z] 



[y] 



Re 



-2 



Figure 4: This contour first goes parallel to the real axis and goes around the cuts at Z + 2m, Z positive 
integer, m — 0,1, . . . 



has no zeroes and its poles are at u^. Being a smooth deformation, H^"^ must also be 
free of zeroes and must have the same number of poles, which are close to the corresponding 
asymptotic positions. Let us denote the positions of the poles by ty^^"^ and w'^^°'\ where 



w 



■w 



m{a) 



(These are indeed far from the real axis and all even cuts.) 
( |5.21| ) can now be written as: 



(5.23) 



u 




(5.24) 



where 



/ p(a)^ 



5.2 Uj related singularities 



x{w 



m{a) \ 



(5.25) 



In the asymptotic solution the singularities of several Y-system elements are given in terms 
of the set {uj}, the set of physical rapidities. In the exact solution, the positions of the 
singular points are smoothly moving away from their asymptotic values and it is possible 
that positions coinciding in the asymptotic limit move away from each other. However, the 
Y-system equations give very strong restrictions also for these positions. 

In the asymptotic solution the set of real zeroes of ^/i"^'^^ coincides with the set of 



physical rapidities for both values of a. After the smooth deformation, the two sets {u^"^}, 
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a = ± may be different, but since are real analytic functions, their zeroes remain 

real. 

Next we use ( |3.21| ) in the Q = 2 case. This gives that Y2 has poles at u'j^^'^ and ^ 

or u^- and ti^-^^ (for all j = 1, 2, . . . A^). Since Y2 is real analytic, both options imply 
that 

uf'^=uf^=u,, j = l,2,...,N. (5.26) 
Now let us use the T-Y relation 



^2^0^2:0 = ^1.07^3,0(1 + 5^2). (5.27) 

We conclude that T2fi has poles at u^~^ and Uj . This is very different from the asymptotic 
solution where q'' = 1. These poles propagate also for higher Tafl functions: (for a > 2) 
the first cuts occur at ±±±11 but there are poles already at u, . 



From 

Y2,al 72,a2?20 



we see that 1"/:"^ has zeroes at uf^ and u"^ . This is a new feature of the exact solution 

l\vw J J 



since T2fi is absent from the denominator of the asymptotic analog of ( 5.28 ). 

Using this last piece of information in the Y-system equation ( ^.181) we conclude that 
since the left hand side has double zeroes at u^, both Y2 and Y^^ must have poles at u^. 

Our final conclusion here is that 

e-^^"^ = 1^ (5.29) 

has zeroes at u^, hence t/;^*-"'* = u'j' and w™"^"^ = uj exactly and we arrive at the final 
result 



Its logarithmic form is 



exp<^ -^Lg^i^gJ . (5.30) 



J^"^=j+jred, (5.31) 

where 

J = In = X; fin - In (5.32) 



'-m 



k=l 



and 



jred = - ^ i^Q * KQy. (5.33) 

(3=1 
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5.3 Simplifying the equations for and 



We can now substitute the results ( 5.32| ) and ( 5.33| ) into the TBA equation ( ^.31 ) and write 

★ s - Li * s + if + j,'ed) i si \ . (5.34) 



1 + y/"^ ' 

11 vw 



(a) 



We will simplify the Y^'^ TBA equations with the help of the following two kernel 
identities. 



-2\j\dvsi{v - u)KQy{w,v) + KQy{w,u) - 25q^i s{u - w) 



Kq{w-u) + 2 dvs{u-v)K^J^{w,v) 



(5.35) 



and 



2 ^ dv si{v — u)f (v) = j{u) 



dv s{u — f ) In 



RmRm I 

Using the above identities and combining ( ^.3C1| ) with ( 5.34| ) we get 



iv). (5.36) 



Yl^'Y^' = Ui M exp<^ 21n 





l\vw 


|21n 


1 + y/,"^ 

l\w 





^s+Y,Lq^[-Kq + 2K^^^s\ 
Q=l 



In 



Rp Rp \ 



RmRm I 

6. Dispersion relation for A 



(5.37) 



* s 



In this section we determine the A discontinuity using the relations (2.9) and the cor- 
responding dispersion relation. At the same time we complete the construction of the 
T-system elements we started in section 4. Our aim is to construct T-system elements 
that are smooth deformations (similarly to the Y-system elements) of the corresponding 
asymptotic variables. As we have seen in section 4, this is possible only partially. If, for 
example, Ta^i is close to the corresponding asymptotic T^f^ solution for the right hand side 
of the diagram, then, in general, even for the left-right symmetric s[(2) cases, Ta-i can be 
very different from the corresponding asymptotic T-system elements on the left hand side 
of the diagram. We will work in this gauge which we call the R-gauge and when we want 
to emphasize this asymmetry we also use the notation T^^ for Ta^s- Of course, there also 
exists an analogous L-gauge, with corresponding T-system elements T^g, which are close 
to the asymptotic solution on the left hand side. 
We start the calculation by writing 



Yi = yi,o 



Ti,iri,_i 

?2,0 



1,1 1,-1 ^^^..^ 



f2,0 



(6.1) 
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The new object here is 



1 



T2,o rf_i Ta,, 



(6.2) 



TjT- is the gauge transformation connecting the two gauges defined above. 



where ta^i 

Corresponding to the three factors in the last expression in ( |6.lD we write 



In- 



where 



A 



R 



In 



A, 



In 



A, 



In- 



hi- 



t 



1,-1 



(6.3) 



(6.4) 



We now rewrite the relation ( |2.9D , using also the results of the previous section in the form 

[A]±,. = Z)Str+Z)[r^±2(,W+£]), (6.5) 

where 



^(±2iV) _ 



and 



A 



{+) 



N 



A 



(a) 



m=l 



In 1 



±{2Af-m) 



D 



{±2N) 
12L 



A 



{-) 



±2N 



N 



m=l 



(a) 



A(^) = In 1 + 



Y 



(a) 



±{2N-m) 

(6.6) 
(6.7) 



So far (except the construction of the T-system elements in section 4) our consid- 
erations are valid for any state of the model. From now on, for simplicity, we restrict 
our attention to states in the s[(2) (sub-)sector defined in appendix C. We think that our 
methods can straightforwardly be generalized to generic states, but some of the subsequent 
formulae are considerably more complicated in the general case. Since the two sides are 
identical in this special case, in the rest of the paper we can omit the upper index of 
the Y-functions Yi"^ , F^"^ or F^"^ . 

± ' m\vw m\w 



In the 5[(2) (sub-)sector in the asymptotic limit we have 



and 



[Af^]±2N = [Ad]±2N = ±2j^^i, 



A 



(0) 



A. 



(6.8) 



(6.9) 



as discussed in appendix B. 

The main part of the calculation is to determine (omitting the upper index R of the 
T-functions) 

r+ 

A,j = ln-^, (6.10) 
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where = and similarly we define Q~ = (T-f^)*. Note that A^(n) = p{u), the 

discontinuity which was left undetermined in section 4. 

Considerable simplification occurs if we express in terms of the T-system ele- 

ments. For the upper sign we find 



■^12^^ = [^^^t^hN + [lnriv+i,i]Ar_i - [lnrAr+2,l]Af- 



(6.11) 



In the gauge we are using Ta^i are meromorphic functions in the strip (—a, a), which means 
that this simplifies to 



D 



{27V) 



.i2fl - [lnr+j2iv. (6.12) 

Next we further specify our gauge: we require that J-~^ has no discontinuities in the upper 
half plane. This is satisfied by the asymptotic function J-^^^^ . We then have 



D 



(27V) 
12R 



In- 



1,1 



[A 



i?j2Af- 



(6.13) 



2iV 



This relation above is in the sense of a discontinuity relation. In this sense one has for any 
function / and q 

f 



In- 



[In/]; 



(6.14) 



provided q is meromorphic around the cut Z. However, as explained in section 2, in the 
sense of dispersion relation we have for example 



dv\hi f]z{v) 



dwfln f{v + — + ie) - In' fiv + - ie)] (6.15) 

9 9 



and the result will be different if we use f /q instead of / when the meromorphic function q 
has zeroes/poles on the cut. The correct interpretation of the original relation ( p.9D can be 
found out from the asymptotic limit of the relation and we find that this subtlety is relevant 
for the = 1 case only and in that relation we have to make the substitution Ai =^ — , 

P2 

where p2 is a polynomial whose roots are the real zeroes of ?2^i. Actually, an equivalent 
(and simpler) way of taking into account this effect is to simply omit the contribution of 
these poles from the dispersion relation. 

Repeating the same calculation for the lower sign we find 



D 



(-2N) 
12R 



[Ing- + ln^+ - lnT+J_27V - [InTjv+i,!] 



+l,lJ-{Ar_l) 



[lnTN+2,i]-N- (6.16) 



Note that the first term is obtained from 

1 



1 



1 

k7 



1 



1 + (n,i)* 



m+i)*(Tri)* T+g- 



2,1)* 



t2,l 



(6.17) 



Adding to the list of gauge fixing conditions the requirement that has no discontinuities 
in the lower half plane (satisfied by the asymptotic Q^^^^) and using the properties of the 
Ta 1 functions we have for both signs 



D 



{±2N) 
12R 



In 



T+ 



[A 



R\±2N- 



(6.18) 



±27V 
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(With the understanding that the correct way of transforming the discontinuity relations 
into dispersion relations means that we have to omit the contributions of real zeroes in the 
= 1 case.) 

At this point we see that the discontinuity relations ( |2.9D (rewritten in the form ( |6.5| )) 
are satisfied by the asymptotic solution. Indeed, all the properties of the T-functions that 
we used in the derivation of (|6.18) are equally valid for the asymptotic T-functions and 



(6.18), combined with ( |6.9| ) and taking into account the asymptotic relations (2.15) and 
the definition of j in ( ^.321) , prove the statement. 

Since in the s[(2) (sub-)sector the two sides of the problem are identical (in the language 
of Y-functions) we have 

[A]±2^ = 2Z)g7)±2(jW+i;i) (6.19) 
and using ( 6.1^ ) this simplifies to 

[A3]±2iv = ±2(jW+j;^jj. (6.20) 
We now define and find the relations 

Yd = ydZd, A3 = Ad + A,ed, A,ed = ln— (6.21) 



Since asymptotically Z^^"* = 1, aI:^^ = we conclude that in the exact solution e^ 



red 



has no zeroes/poles (and has constant asymptotics for large |u|). This means that if we 
define (analogously to F^") (n) of section 5) 

T,,^{u)= i dvi^',,^{v)K{v,u) (6.22) 

and write the corresponding dispersion relation we find 

oo 

PAred = 2 jiSd * iK^'""^ - ^'-'^^ (6-23) 

N=l 

Here the left hand side is the contribution of the narrow strip only (the pole terms vanish) 

and the right hand side is simply the sum of contributions of the integrals along the even 

cuts. Using the method explained in section 5 we can "integrate" this relation and get 

oo 

Ared = 2 J2 & * (^''^' - (6-24) 

N=l 

6.1 Determination of A/j 



The dispersion relation for Ar can be written down using (|6.6|) and (|6.18D . Proceeding as 
before we have 

r^(n) = / dvA'ji{v)K{v,u) = {VAr){u) + <f dv [In T+^iv) - In {v)]K{v,u) 

Jto i[7,2,oo] 

+ / d^; [ln'T+-^{v) - In {v)]K (v , u) . 

(6.25) 
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The alternative ealculation uses the expUcit form of D^"^^ and its direct integration along 
the cuts. The contribution of the upper cuts is 
oo „ 

N=l 

= - j dvK^'^\v,u)A^2^'{v)+ (f dvK{v,u)A'_{v) (6.26) 

J — OO J [7,2,00] 

00 /"OO ^ p 

™_i J —00 . 1 JfT.m.ool 



Let us consider now the pole terms {§ terms) separately: 

00 

Tl -'[7+l,"''+l-)0o. 



m=l "'[7+l,"»+l,oo] 
00 „ 

+ J2f dv ^[-+11 {v, u) In' {v) 

00 „ 

-J2f dvK^^+'\v,u)\n'Tm+i,i{v) 

00 „ 

J2 f dvK^^-'\v,u)\n'Tm+i,i{v) 



(6.27) 



Wc first note that in the integration contours m — 1 can be changed to m + 1. This is 
possible because Tm,i has no zeroes/poles/discontinuities on the m — 2 cut for m = 3, 4, . . . 
and T2^i has no poles/discontinuities on the real cut while the contribution of its real zeroes, 
as we have seen, has to be omitted. We can also change all integrals to have contours of 
the type [7, m + 1, 00] but here the change of 7 it 1 to 7 requires to add to this sum the 
corrections 

00 „ 00 „ 

- V i dvK^^-'^{v,u)WTm+i,i{v) + V i dvK^^+^\v,u)WTm+i,i{v). 

777=1 •^[7,7+1] „=1 ^[7-1,7] 

(6.28) 

The rest of the sum (after taking into account that most terms cancel) simplifies to: 

* dv K{v,u) In' T2,i{v) - f dv K {v , u) In' T^^{v) + f dv K+{v,u) In' Ti^i{v). 

(6.29) 

Adding also the pole terms corresponding to A_ 

/ dv K{v, u) [in' (v) + In' Tf ^ {v) - In' T2 1 {v)] (6.30) 

Jb,2,oo] L ' ' 'J 

the sum of all pole terms becomes: 



00 „ 

- V i dvKl^-'\v,u)ln'Tm+i,i{v) 

^=1 ^[7,7+1] 

+ V/ dvK^"'+^\v,u)ln'Tm+iiiv) + (f dv K{v,u)\n' T+^iv). 

m=0 •^[7-1,7] ' -'[7,2,00] 
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This can be simplified furtlicr using the relation Tm,i = {—l)"^~^Trn(^rn and the facts that 
am has no zeroes/poles near the physical strip and has no zeroes/poles on the even 
cuts: 



oo „ oo „ 

+ f dvK^''\v,u)ln'T[^^^\v)+ (f dvK{v,u)ln'ativ). 

J—oo J [7,2,00! 



(6.32) 



It remains to calculate the parts. We write these contributions using the combinations 
(already used in section 4 for the construction of the T-system elements) 



=Sf± = ln 



-T2 1 



Y± 



TmTm+2 1 + 



m\vw 



(6.33) 



and find 



■OO m=l -^-^ 

/OO °o /•OO 

-OO r„=l -^-oo 

/.OO °o „oo |- 

+ / d?; W (t;, u) In' r^'^J (v) + / di; ^1"^+^] (?;, u) In' t|^1 (?;) + In' r|]'^2 



m=l 



(6.34) 



Adding everything together and using the identity r^r^ = 1 in the form In' Tm + 
In' Tm'^^^ = we finally have 



r^PP«'"(u) = / dv K{v, u) In' a+{v) 

J[7,2,oo] 

Similarly for the contribution of the lower half plane cuts we have 

ri^^e'^(n) = / dv K{v, u) [W a+{v) - In' jr+ (v) - In' (v)] 

i[-oo,-2,-7] 

/OO °o „oo 

dvK^-'\v,u)^l~'^\v)-J2 / di;i^[-"*l(?;,u)^4(?;). 
-OO ™_i J—oo 



(6.35) 



(6.36) 



6.2 The XXX gauge 

To calculate /S.r we need to fix the gauge completely. We will call this complete gauge 
fixing the XXX gauge and it will not be our final gauge choice yet. Later we will have 
to perform an additional gauge transformation to arrive at the gauge in which all the T- 
functions are smooth deformation of the asymptotic solution. This "intermediate" gauge, 
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which win be used in this subsection, we call the XXX gauge because in the g — )• limit 
the system of functions Y^^^^^, Tm,i decouples from the rest and they form the Y-system 
and T-system of the XXX spin model (the latter in a gauge most natural in that model, 
namely in which the Tm^i functions are polynomials in the spectral parameter u). 

The XXX gauge is fixed by adding the requirements (in addition to those already 
imposed in section 4 and this section): 

• ^xxx- zeroes/poles in the upper half plane (no zeroes/poles on the real line) 
except poles at positions (j = 1, 2, . . . , N). 

• Gxxx'- zeroes/poles in the lower half plane. 

Recalling that T^^^ = TicJi we can express VAr comparing ( |6.25D to the sum of 
dOSD and ( lO^ ). We find 

J —oo J —oo 

oo „oo „ 

dv^;^{v)km{v,u)+ f dv [W T+iv) -In' T+{v)]Kiv,u) (6.37) 



+ f dvK{v,u)[-\n'T+{v)]. 

J[-oo,-2,--y] 

We can calculate the pole terms (^ terms) of this expression using the residue theorem: 

N N oo ^ . . A ■ \ 

-2mJ2K{u+,u) + 27TiY^ Yl \K{uj + —,u) - K{u- + —,u) \ . (6.38) 

j=l j=l v=-oo yd 9 ) 

Their contribution, after "integration", as explained in section 5, becomes 

N N oo ^ . . A ' ^ 

j=l j=l v=-oo ^ ^ ^ ' 

Note that using the definition 

^tiu) = n n si .^+^ ' ^'-''^ 

j = lu=-oo "^l^'^i + g ) 

where the notation indicates that this function has exactly the same zeroes/poles as , 
we have 

N 

e^i^) =^+(u)YlS{u,up. (6.41) 

i=i 

We can now write the "integrated" version of (|6.371 ): 

Aniu) = - r dv^["'\v)K^^\v,u) - r dv^]~^\v)K^~'\v,u) 

J — OO J — oo 

(6.42) 



^ / dv ^^{v)k^{v,u) +yV{u) -In x'^{u). 
m=i 



-31- 



Note that this formula is valid just above the real line in the strip < Imn < ^. 

The last term in (|^) is necessary to balance the large \u\ asymptotics of the equation. 
We know that in this limit (just above the real line) 

.i4(u), ^_(n) w -In |n|, A^(u)«ln|n| (6.43) 

and we see that (l6.42| ) is asymptotically correct because if a function / behaves asymptot- 
ically as f{u) w In then convolution with the kernel K^"^ gives 

/^i^H ~ _iln|n|. (6.44) 

Finally by manipulating the integrals containing .if_ we write our final result for 
the discontinuity Ar: 



Ar{u) =^-{u)- j dv[^-{v) + ^+{v)]K{v,u) 

^ f'OO 

dv ^m{v)km{v,u) +W{u) -In x'^{u). 



(6.45) 



Note that the last term is chosen here so that Aj{(ib2) = is satisfied. The price we have 
to pay is that this term creates a cut along the imaginary axis (but neither A^(u) nor 
Q^ni'^) have discontinuities there). 

Our final result for the discontinuity A(u) in the s[(2) (sub-)sector is 

A = 2Ar + Ad + Ared 

°° -r)+ (6 46) 

= - 2(if_ + ^+) i - 2 ^ ★ + 2W - LlnxV In + A,cd- 

m=i (^1 )* 

Here we have used ( C.30| ) to write 



Ad = -Jlnx2 + ln-^. (6.47) 

v^i )* 

We note that the renormalization J— )-L = J + 2is due to the addition of the zero mode 
contribution in ( |6.42 ) and it is universal (we see this here for the case of states in the 
s[(2) (sub-)sector). Indeed, this universality was shown for generic states in The 
physical meaning of L is that it is the maximal value of the J-charge within the psu{2, 2|4) 
supermultiplet to which the given state belongs. 

The result ( |6.46| ), together with calculated in section 5, completes the set of TEA 
integral equations of section 3 for the excited states in this sector of the model. This system 
of integral equations is now closed. It still has to be supplemented by the quantization 
conditions for the discrete parameters appearing in the source terms. This will be discussed 
in section 9. 

At this point we would like to emphasize that the Y-system and discontinuity relations 
(with some additional qualitative information on local singularities) determine the set of 
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TBA equations completely. We can summarize the logic of calculating the full A (including 
the dressing phase part A^) as follows. First we fix a gauge such that the T-factors in ( |6.1| ) 
satisfy ( |6.18| ) and thus the discontinuity relations simplify to ( 6.20| ). Next we take ( |6.21| ) 
as an Ansatz and use the (nontrivial) computations presented in appendix B to show that 
(6.20) is now further reduced to 



[A, 



red\±2N 



±2i 



red^ 



(6.48) 



which can be solved easily using ( 6.22 - 6.24 ). Of course, by writing ( 6.21 ) we actually use 



the known formulae for the dressing phase part and verify it satisfies the relations ( B.23 - 



B.24). This simplifies our job here. However, even if we had not known the solution for 



Arf given by ( |6.47D we could have calculated it from the discontinuity relation (|B.23D by 
transforming it into a dispersion relation as we did for the other building blocks. 

In this section we obtained the result ( 6.45| ) by a long direct calculation. An alternative 
logic could have been to simply postulate the result ( 6.45| ) for p{u) = Ajj(n) and show that 
the requirements on the upper/lower half plane behavior of the jQ^ functions, which 
we used in the calculation, are indeed satisfied. 

To show this, we start from 



0"! 



-a 



(0) 



(6.49) 



where 



af = exp \ ^1^1 * {k- -K+) + ^^m* {k 



[-m-l] _ ^[m+l]^ 



5.50) 



m=l 



and q is defined in ( [4.22| ). Here we introduced the notation 



k[v, u) = : 



27ri V — u 



-a\ ,.Ja 



.51) 



Note that K{v,u) — k{v,u) is regular at v = u. 
The discontinuity relation is 



q{u + ie) 
q{u — ie) 



-p{u+ie) 



\u\ > 2 



(6.52) 



and just above the real line, in the strip < Imu < ^ we have 



p = - (if_ + ^+) ★ K - ^ ^„ * + W - In; 



(6.53) 



and 



m=l 



.(0)- 



exp J + ^1'^ ★ {k— -K)+^^rn* i^^' 



-m-2] _ ^M) 



.54) 



771=1 
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and finally 

-^^xx = {T^n''--' = $ exp {(^_ +^^)iK 

(6.55) 

m=l 

From this result we see that there are indeed no zeroes/poles/discontinuities in the upper 
half plane, except for the poles at the positions . 

Similarly, just below the real line, in the strip — ^ < Imu < we have 



= exp <j ^_ + ^1^1 * (k - + ^ ★ - ) j> , (6.56) 

m=l 



p = - (^_ + ^+) i - ^ ^„ * A;^ + W - Inx^ (6.57) 

m=l 

and noting that from ( 2.11| ) we have / = — ^+ 

oo 

TO = / - p = + (^_ + ^+) i iv: + ^ ^„ * A;„ - W + lna;2, (6.58) 

m=l 

and finally 

^ixx = {T^n'^-"' = 7^ -r { - + -^+) ^ ^ 

oo 

* (k - K++) +Y,^m* («[""! - - - k[™+2])|. 



.59) 



m=l 



We see that there are indeed no zeroes/poles/discontinuities in this expression in the lower 
half plane. 

6.3 Complete gauge fixing 

It is possible to show that the XXX gauge we have been using in this section is a complete 
gauge fixing. This means that given an exact solution of the Y-system equations, which 
also satisfies the reality conditions and the relation between y_ and 1+ and has the cut 
structure described in section 2, we can always construct the corresponding XXX gauge T- 
system elements T^^^, which satisfy all the requirements listed below and this construction 
is unique. 

The requirements are 

. T^XXX _ 1 

• -^0,s = ^■ 

• T^^^ are given as constructed explicitly in section 4 using the chain lemma. 

• T^f^^ have discontinuities along the cuts =b(a + 2m), m = 0, 1, . . . 
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• The roots of T^^^ in the physical strip are the same as those of and there are no 
poles in this strip. 

• ^xxx- ^° zeroes/poles/discontinuities in the upper half plane (no zeroes/poles on 
the real line) except poles at positions (j = 1, 2, . . . , A^). 

• ^xxx- zeroes/poles/discontinuities in the lower half plane. 

• The large \u\ behavior is the same as for the asymptotic solution. 

We now perform the final gauge transformation which brings the T-system solution 
in a gauge (we will call it the BA gauge) where the T-system functions Ta^s for s > are 
close to the asymptotic (Bethe Ansatz) solution ri^i . This is achieved by defining 

We see that Ti^i has no zeroes/poles in the physical strip. From 

-^^ = ^ -^^xx, = V ^^^^ ^^-^^^ 

we see that IQ~ has no zeroes/poles/discontinuities in the upper/lower half plane. 
From the list of requirements above and the modification induced by (6.60) we see 



that the BA gauge T-functions are indeed smooth deformations of the asymptotic solution, 
satisfying the same requirements. Note that the final results for p = A/j and the full A are 
unchanged since the transformation ( |6.60 ) is meromorphic. 



7. Simplified Yi equation 



In this section we want to simplify the TBA equation ( 3.30| ) in order to be able to compare 



it with the results of [p4|. First we simplify Ar and Ared and use these results in the Yi 



TBA equation (3.30). 



7.1 Simplifying 

Let us introduce the new TBA variables Wa with the definitions 

Ya = T^Wa, a = 2,3,..., Y^ = Wi. (7.1) 

Wa {a = 1,2, ... ) have no zeroes/poles in the physical strip. In terms of these variables 
we rewrite ( [3.29D in the form 

lnWa = {2^a-i + lnWa+i + lnWa-i-La+i-La-i}*s, = 3,4,... (7.2) 
In W2 = {2ifi + In W3 + In VFi - L3 - Li - In rf } * s. (7.3) 

The simplification of ( |6.45|) is based on the kernel identity [^6| 

ka+5a,lS = S-k{ka+l+ka-l), a = 1, 2, . . . (7.4) 



-35- 



where ko{w,v) = 20(4 — w'^)K{w,v). The convolution of (7^) with {ka + ka-2) (from the 
right) gives for a = 3, 4 ... , after using ( |7.4D 

In Wa * (ka + ka-2) = (2^a-l + In Wa+1 + In Wa-1 - La+l - La-l) * ka-l (7.5) 

and from ( ^^ ) we get 

In W2 *{k2 + ko) = In W2 + (2^i + In W3 + In VFi - L3 - Li) * A;i - In rf * ki . (7.6) 



Using the methods we employed in section 6 for the calculation of pole terms we can 
evaluate the convolution In * /ci : 

N 

Inrf + fci = 2>V->Vo, >Vo(n) = J^[a(tx, u+) + a(n, nj)]. (7.7) 

i=i 

The equations (^) and ( [7^ ) can be used to express the terms containing hi ( 6.45| ). 
Many terms cancel and we find 



2^^rn*km = ^iLQ+ LQ+2) -k kq - hi W2 
m=l Q=l 

+ 21nTy2 * A' - InFi ★ fei + 2W - Wo. 
This can be used to write the simplified Aji formula 

2An = - f; (Lq + Lq+2) * A:q + In (1 - ^) | 



(7.8) 



1 \ / 1 



(7.9) 



21n<;y2( 1-^ 1 ( 1-^ ) ★i^ + Wo-21nx2 + lnyi*A:i. 



7.2 Simplifying Arcd 
We start from 



00 „ 

Ared(n) = 2 [hdvj,,^{v + ie)[K^^''\v,u)-K^-'''\v,u)] (7.10) 

^ /"OO 

= ^ / dw Lq{w)J:q{w,u), (7.11) 

^ 1 J —oo 



where 



oo „ 

^q{w,u) = 2^2 fdv [k^^\w,v + ie) - K^-^\w,v + ie)^ i^K^^^\v,u) - K^'^^\v,u)^ . 

N=l •' 

(7.12) 

The kernel identity we need here is 

(7.13) 
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Using it we find 

^ ( (-2 



Y.q{w,u) = 2^2 \ / dv KQy{w,v)k2Niv,u) - k2N+Q{w,l 
N=l 1-^-2 

= 2Iq{w,u) -2 J dvKQy{w,v)Io(v,u) = 2kQ(w,u). 



(7.14) 



7.3 The Yi equation 

The results of the previous two subsections can be used to simphfy the full discontinuity 
A. We get (with L = J + 2) 



A = ln^ 
(P: 



|_ + A(^) + A(i2) + y^^ + In jy^ (^1 _ I _ Llnx^, (7.15) 



where 

oo oo 

A(^) = Y,Lq^ {2K^ -kQ)-Y,LQ* kQ.2 (7.16) 
Q=l Q=3 

and 

A(12) = inYiicki - 21n j^s " j^) " } * ^- (^-1^) 
The simplified TBA equation for Yi becomes 

InFi = Injys (^1 - |*s-L2*s-Ais (7.18) 



= lnjy2(^l-|^J I *s-L2*s-Aeff (7.19) 
where 

Aeff = A^-^) + A(^2) ^ ^(source) _ ^ 1^ (7.20) 

and 

7-)+ 

^(source) = >v + In ^1 . (7.21) 
Using the results of appendix B we can write 



= iW'^' [^ho(-.4) - G(>'.-7)l} . ("2) 



where 

iG{u,0 = '^(^-^,^^-<^{xiu),0+i^iu,0 (7.23) 
and the source term simplifies to 

^(source) ^l^J^I^^ +2X;{G(^,4)-G(u,a;7)}. (7.24) 
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If we use for G(n,^) the alternative representation 



Giu, = r dv In [ ] Io{v, u) (7.25) 



2 



x{v) — ^ 



we can check that our simphfied equation ( [7. 191 ) agrees with eq. (4.10) of ref. p^ . 
8. Canonical (and hybrid) equations for Yq 

Using the Wa variables introduced in section 7 we can rewrite the Y-system equations 
(|3^2lH3^ ) in the form 



° " = Ca, a = l,2,..., Wo^l, (8.1) 

where 

In = 2if„_i - La+i - La-1 - 6a,2 In a = 2, 3, ... , (8.2) 

In^i = 2if_ -La. (8.3) 

Using the chain lemma of section 4 we can transform (|8.lD into the integral equations 

oo 

lnWa = Y,lnU*ia-^*Ka. (8.4) 

A=l 

Here A is given by ( |6.46|) . 
8.1 Kernel identities 

We now list a number of kernel identities that are needed to simplify the integral equations 



(|8J). 



We start by writing the kernel defined by (|A.18|) as 



^) = 1^Kq+m{u -v)- ]-Kq_m{u -v) + ^ KM-Q+2j{u - v) 

^ ^ i=i (8.5) 

+ lA:Q(n,.M)-iA:Q(n,.[-^^^]). 

Next we write the identity 

K{y, v)Ka{v -u) = K{y, u^^^)K{v, ^x^) - K{y, u^-^^)K{v, txl""]) - K^iy - u)K{v, y) (8.6) 
and integrate with respect to v just above the cuts using the result 

j\dvK{v + ie,a) = (8.7) 



and obtain 



KiKa = -Kya - -Ka. 



-38- 



This can be used to get the further identities 
and for the "fermionic" kernels defined by ( A.13| ) 



Kl^ = -Kya + -Ka = Ka + K i Ka, (8.10) 
= ^Kya -^Ka = Ki Ka. (8.11) 



The most important identity is 
which is proven in [^6|. Using also 

iQ+i+(Q~i = KQa (8.13) 

we can write 

^Q+i ^ + ^ ^ + 2K^ iKa = KQa + 2K^a = -K%y (8.14) 

Finally we write the cliciin leninici for ha — -^r^. Since it satisfies 6^ = fta+i^a-i hi this 
case = 1 and the only nontrivial object is 

h+ 

p = ln-^ = -lnx2 (8.15) 
(^1 )* 

and we find from the chain lemma 

4 = -ln6a = -Inx^ ^ (8.16) 
8.2 Canonical TBA equations 

Collecting all the terms proportional to or Lq and using the above identities and 

the simplifications of subsection 7.2 we can rewrite ( ^.41) . This can be called the canonical 
Wa TBA equation and is of the form 

CXD OO 

In Wa = -L£a + /a + 2 J] ^m*KZl + 2if- * Kl^ + 2Jf+ iKf + Y, *^'?(2) • (8-17) 

m=l Q=l 

Here an alternative form of the .J^± terms is 

(8.18) 
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and the source term /„ can be written as 

fa = - Inrf * {Ka+i + Ka-i) - (2W + In j i Ka- (8.19) 

The source term can be simpHfied using the chain lemma for 

N N _^ 

-Daiu) = n St^,^{u,U^) = n {^i2)(%,«)}" • (8-20) 
j=l 3=1 

T>a satisfy 



V+V- = Va+iVa-i, a = 1,2,..., Vo = \. (8.21) 

However, D2 has poles at and zeroes at uj , so the chain lemma does not directly apply 
for the T>a functions. This problem is solved by introducing 

Q[-+^(.)Q[-il(n) _^ 1 

«aW g[-a-l](^)Q[l-a](^) ii S,^^{u - U^)Sa-l{u - Uj) ^^'^^^ 

and noting that = ^ also satisfy the chain lemma equations with = 1 and do not 
have any singularities near the physical strip and thus can be represented as 

P„ = exp|-ln^^iK,|. (8.23) 

Using this representation the source term becomes 

N 

fa{u) = - 1; In Si*^^) {u^,u) - [In rf * (i^„+i + i^„_i)] {u) 

N 

+ [In Sa+l {U -Uj)+ In Sa-1 {U - Uj)] - 2[T' ki i Ka] {u). 



(8.24) 



Here we used the result 

W = T'i^ki, (8.25) 
which can be proved by using the residue theorem. We have defined 

N 

T^%u) = Y\nt(u-Uj±ie). (8.26) 

It is easy to prove that 

TV 

[{T' - T-')i.e^ {u) = -Y, [ln5„+i(tx - Uj) + \nSa-i{u - Uj)] , (8.27) 

i=i 

{T' + T-')i.ll = In rf il (8.28) 
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and with the help of these relations and the results in the previous subsection we can write 
the final form of the source terms: 



N 

faiu) = -^ln5i^»)(tx„n) 

N 

= -^ln5i^»)(u„n) -2 [T^^i^i^,] (n). 



(8.29) 
(8.30) 



8.3 Hybrid TBA equations 

We can get rid of the infinite sums containing the functions by the following trick 
We write the ( |3.25 ) TBA equation as 



= {rm+i + rm~i " L^+i) *s + 6^,1 In 



1 -y_ 



1-n 



*s, m = l,2,..., (8.31) 



where WG introduceci the notation Vfji = ln[l+ym|„^], m = 1, 2, . . . , ro = 0. We now assume 
that the kernel functions J^m satisfy the relations 



J^m — S-k (J^m+l + '^m-l) — jm, m — 1, 2, 



^ = {i 

with some jm- We take the convolution of ( |8.31| ) with and sum over m. We find 



^ ^ •^m * — ^ ^ L/m, 'k S rj^fm—l ~\~ ^ ^ * Jm ln 



r7i=l 



m=l 



m=l 



We now choose = K^^. In this case 

and using the above trick the TBA equations can be brought to the form 



.32) 



.33) 



.34) 



In Wa = - L£a + fa + '^Ta-l -k s + 2ri -k s i K 



ya 



2 In 



1 -y_ 



i s k KZ, + 2if_ i + 2^+ ★ KT 



.35) 



Here K^!^^ = is understood. Our final result is in complete agreement with the corre- 
sponding results valid in the special cases studied in p3] and 



9. Quantization conditions and exact Bethe-Yang equations 

In this section we formulate the quantization conditions and the exact Bethe equations 
which determine the discreet parameters £,m,j, £,m,j and uj occurring in the source terms of 
the TBA integral equations. 
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9.1 Quantization conditions 

In this subsection we quantize tlie roots occurring in the functions 

Mm J^m 
Tm{u) = Wt{u- Cm,j), f^iu) = ^(^i - Imj), 171 = 2,3,... (9.1) 

j=l j=l 

(Note that fi(n) = 1 and the physical rapidities uj occurring in Ti(n) = n*('" ~ ^^^^ 
be quantized by the exact Bethe equations discussed in the next subsection.) We assume 
that (as is the case in the asymptotic solution) the functions ^^'^^ poles around 

the points ^mj or Cm+2,j and therefore the zeroes on the left hand side of the Y-system 
equations ( ^3.17 ) must be accompanied by corresponding zeroes also on the right hand side. 



This leads to the quantization conditions 

l + =°' rn=l,2,..., j = l,...,Mn^+i. (9.2 

These are well-defined even for m = 1 since 

(^2,i + ie) 1 _ y_ (^2_^. + ie) 



(9.3) 



f±2l 

Similarly we assume that the functions Y^^J have no poles around the points or 
Cm+2,j and therefore we have the quantization conditions 

l + Y^^JU+ij) = 0, m = l,2,..., j = l,...,Mra+i. (9.4) 

These are well-defined even for m = 1. 

We discuss an important special case in detail. In this special case (which is relevant 
for example for the case of twist-two states in the s[(2) sector) we have no roots at 
all and all roots are real. Moreover all Afm are even numbers. Then we can make the 
following definitions. The (|3. 25 - 3^61) TBA equations can be written 



Ym\vw = TmTm+2 exp {fra * s} , (9.5) 

where 

, ^ [l + ^m-l|^,^(^)][l + ^m+l|^,^o(^)] ^ , 2,. 1 f ^-Y-{u) \ 

/ml") = In T—TP r-. + 9(4 - u )6m,i In -— - (9.6 



and we define 



Bmin) = f P r dv ■^7^''\ . (9.7) 
4 iu-v)gn 



oo sinh ■ 2 



Note that the principal value prescription V is actually superfluous at the points we will 
need this function since 

/m(Wij) = 0, m = l,2,..., j = l,...,Mm+i. (9.8) 
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We also define 

R^{u) = 2 ^ arctan tanh ~ ^^'^'^^ . (9.9) 
i=i 

In the special case the quantization conditions take the form 

Rm{im+l,k) + Rm+2{im+l,k) + ^m{im+l,k) = '^'^^m+l,k m = 1, 2, . . . , /c = 1, . . . ,Mm+l, 

(9.10) 

where 

mteger it is oad, 

(9.11) 



^m+l,k 



halt — integer it is even. 

^ 2 



We note that the above quantum numbers fm+i^k are not free parameters. They are 
determined by the state under study and in principle can be calculated by computing the 
left hand side of ( 9.10| ) in the asymptotic limit. 

9.2 Exact Bethe equations 

The quantization conditions for the physical rapidities Uk are the exact Bethe equations: 

= -1, A: = 1,2,...,A^, (9.12) 

where the analytic continuation (denoted by the subscript (*)) means that we have to 
analytically continue the function from the real line just below the cut line at — | coming 
down between the branch points ±2 — |, and then going back to the real axis through the 
cut. Using our previously introduced notation, for any function F[u) we have 

F(,)(n) = ((F-)J+(n-ie). (9.13) 

Our definition is in agreement with the transformation rules {u) = xf (u) and 

'^si'(2^){*) (^i' ^) ~ '^s['(2Y ' ^-^ ■ following transformation rules 
hold: 

f = FiKi f^^^^u) = f{u)-F+{u-ie), (9.14) 

i, = FiKi ij^,){u) = 'il^{u) + F+{u-ie). (9.15) 

Using the last formula we can write the following further identities 

F i Kl') = F i Kl\^^, (9.16) 

FiRyj) =Fi Kl). + i^i'"'^ (9.17) 
/ (*) ^ 

(F i Ky,)^^^ = Fi + (9.18) 

Using this set of identities we can rewrite the exact Bethe equations: 

In Yi(^)(ufc) = 27rii/fc, i^k ■ half - integer, k = l,...,N, (9.19) 
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where 



InY, 



L ln(x+x ) + /i(^) + 2ri * s^' + 2ri * s ★ 

1 - y: 



2 In 



oo 



(9.20) 



and 



N 



We note that since 



2^2 



-2iR2 + In 1 



1 



(9.21) 



(9.22) 



and here the second term vanishes at u^, the first term in the third line in ( p.20 ) can be 
substituted by — 2zi?2- We also note that the e prescription in the third term in ( |9.20| ) is 
not really needed since ri vanishes at u = Uk- 

In ( |9.19| ) i^k are momentum quantum numbers and can be used to label the state 
(instead of the particle momenta). 



Finally we give an alternative variant of ( |9.21 ) using the regularization introduced 



m 



^ f°° 0(v) 



N 

2j;in 



(n,-n++) ' "'"^"^ 



Xj — x^[u) 



(9.23) 



At u = where we need it, we have 

N 



/iwK) = - ^ln5iy;(n,-,^fc) + 2 r d^; kH^^^ 



^{v + ie,Uk) 



N 

2j:in 



{Uj - Uk ) 



2i , 



9 X, 



(9.24) 



Again, the final formula for the exact Bethe-Yang equations agrees with the results 



obtained in [34| by contour deformation techniques. 



10. Simplifying the energy formula 

The energy of the multi-particle state we have been studying in this paper is given by the 
formula 



j=l Q=l'^-°° 



(10.1) 
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where 

£{p) = Jl + Ag^sin^?-, e^^ = 4^ (10.2) 

gives the energy of a physical particle with momentum p in the string model and p'^{u) 
is the momentum of the bound state Q-particle with rapidity u in the mirror theory. Not 
only the above energy formula, but many other results in this paper contain an infinite 
sum of convolutions of the form 

oo 

Y^LQ^fQ. (10.3) 

Q=l 

The results that contain an infinite sum of this form include the hybrid equation for Yi, 
the exact Bethe equation and the results for both Y_/Yj^ and Yj^Y_. The fact that allows 
the simplification of the infinite sum is that in all these formulae the coefficient functions 
/q satisfy the functional relation 

/q+/q =/q+i + /q-i- (10.4) 

We find that if we express Lq in the sum in terms of the T-system functions Tq^ then using 
the above identity most of the terms cancel and we end up with a simple finite expression 
containing the single function Ti^ only. (In some cases two functions remain, Ti^ and 

72,0.) 

For the case of the energy formula we have 

fQ{u) = ^P«(n) = ^ (x[«l'(n) - x[-«]'(n)) (10.5) 
and the full energy expression can be written as 

ig J. , , ^ , , u 



where the contour of the u integral has to lie a little above the real line. 

We discuss just one more example here. In the case of the ratio the relevant 

coefficient function is 

/q = K^Q\ - K^-Q] (10.7) 

and the simplified final formula is 

^exp|21nri,o * i^t'lj . (10.^ 
1,0 ' J 



Y— RpBfji 

y+ BpRmT^ 



In our opinion the fact that the energy formula can be rewritten such that it depends 
on the single variable Ti^ only might indicate that there is a transfer matrix formulation 
behind the exact TBA equations in the model. This may also be an important step towards 
the NLIE description of the system since Ti^ can easily be expressed by the elementary 
objects appearing in that approach 54]. 
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11. Conclusions 



In this paper we derived the TBA equations for the 51(2) (sub-) sector of the AdS^ x 
mirror model based on the Y-system and the discontinuity relations proposed in [44|. 
The proposal for the discontinuity relations was based on the analyticity properties of the 
solutions of the ground state TBA equations and was conjectured to be state independently 
valid for the excited states as well. 

In this paper we have shown that the discontinuity relations hold nontrivially for the 
asymptotic solutions of the excited states. This corroborates the state independent nature 
of the discontinuity relations. In addition we studied the discontinuity relations carefully 
and concluded that a technical subtlety requires the use of a refined interpretation when 
translating them to dispersion relations. 

Our derivation of the TBA equations is based on the fact that the Y-system equations, 
the discontinuity relations plus some qualitative information on the local singularities of 
the Y-functions and their asymptotics at infinity together make possible to transform the 
functional equations into TBA integral equations in a unique way. 

In this derivation we assumed that the Y-functions for the excited states are smooth 
deformations of their asymptotic form, so the qualitative information on their local singular- 
ities and on their behavior at infinity can be read off from the explicitly known asymptotic 
solution. 

Since as we have proven the asymptotic solutions satisfy the limiting functional equa- 
tions exactly, by construction they also satisfy (the Yg — t- limit of) the TBA integral 
equations and the quantization conditions, including the Bethe-Yang equations. An im- 
portant consequence of this observation is that the asymptotic limit of the exact Bethe-Yang 
equation ( |9.2[1| ) accounts for the Beisert-Staudacher equations ( p.31| ). This fact has not 
been proven analitically so far, though it was an important starting point in the 5-loop 
tests of the mirror TBA equations |41, 42, E^]. 



Beyond the derivation of the TBA equations we also constructed the T-system elements 
in a special gauge in terms of the Y-functions. The benefit of this construction is 2-fold. On 
the one hand the discontinuity relations are a lot simpler in the language of T-functions 
and their introduction makes the derivation of the TBA equations easier. For example 
with their help we could show that to derive the TBA equations the knowledge of local 
singularities of the Y-functions lying only within certain finite strips around the real axis 
is needed. On the other hand we recognized that important infinite sums of the TBA 
problem simplify drastically if they are expressed in terms of the T-functions. The most 
important such simplification appears in the energy formula, which can be expressed as a 
simple expression of a single T-function Ti^. This might indicate that there exists a hidden 
transfer matrix formulation of the model. 

Independently of these speculations, we think that our construction of the T-system 
and the expression for the energy in terms of Ti q gives an important step towards the 
NLIE formulation |51, 52] of the AdS/CFT spectral problem. There the T-functions are 
more fundamental objects and can be expressed easier than the Y-functions since the NLIE 
construction is based on the T-Q relations of the model. For example in the approach of 
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1 53] and the Ta^s functions are expressed by the Wronskian determinants of certain 



fundamental Q-functions whose combinations are the unknowns of the NLIE. 

As a final remark we note that recently |55] with the help of the T-Q relations the left 
and right SU{2) wings of the TBA equations could be resummed by a hybrid-NLIE and 
so the number of unknown functions were remarkably reduced. 
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A. Notations, kinematical variables, kernels 



In this paper we adopted the definitions and conventions of ref. |34]. For completeness, in 



this appendix we collect these definitions and give a list of all kernel functions used in the 
paper. 

We will use the notation f^{u) = f{u it -) for any function / and in general f^°'\u) = 
f{u + |a). We will also use w'^ = w ^ ^ for w some parameter. 

Most of the kernels and also the asymptotic solution of the Y-system is expressed in 
terms of the function x{u): 

x{u) = -{u — i\l 4 — ti^), Imx(u) < 0, (A.l) 

which maps the u-plane with cuts [— oo, —2] U [2, oo] onto the physical region of the mirror 
theory, and the function Xs{u) 



x,(n) = ^ l^l + yi- -ij , |x,(n)|>l, (A.2) 

which maps the u-plane with the cut [—2, 2] onto the physical region of the string theory. 
Both functions satisfy the identity x(u) + = u and they are related by x{u) = Xs{u), 
and x{u) = l/xs{u) in the lower and upper halves of the complex plane respectively. 

The momentum and the energy £q of a mirror Q-particle are expressed in terms 
of x{u) as follows 

i i ~ x(^u — -Q) 

PQ = gx{u- -Q) - gx{u + -Q) +iQ , £Q = log— 1—-. (A.3) 

9 9 x{u + ^Q) 

Three different types of convolutions appear in the TBA equations. These are: 

/oo f2 
du f{u) )C{u,v) , fi)C{v)= / du f{u) )C{u,v) , 
oo J -2 

fi)C{v) = (/ ' + du f{u) )C{u, v) . (A.4) 
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nel /C depends on a single variable, then the convolutions in ( [A.4p are understood as 
f duf{u)IC{u — v). For a kernel /C and parameter a G M we often use the notation 
IC^''\u,v) =IC{u+f,v). 

The different but equivalent formulations of the mirror TBA equations^^ contain a 
number of kernels which we list below. 

We start with kernels depending on a single variable: 



2 



■s(^i) = TT- 3- log * (^) = 1 ' *(^) = tanh[— u] 

ivr? du 4 cosh 4 



I d ^ ^ , , 1 gQ 



u 



Kqiu) = £ log Sq{u) = ^ ^ , , Sq{u) - ' 



2iTi du IT Q"^ + g'^u'^ ' ^ _|_ ' 

1 d ^^'^ 

Kmn{u) = ^ ^ logS'MAf(ti) = Km+n{u) + Kn-m{u) + 2 ^ KN-M+2j{u) , 

M-1 

Smn{u) = Sm+n{u)Sn-m{u) Y[ SN-M+2j{uf = Snm{u) . (A. 5) 

i=i 

The fundamental building block of kernels which are not of difference type is: 

zvrz du zm -y/4 — u — v x[u)x[v) — 1 

An important function in the equations for J^") and A is cr{u, v) which is the logarithm 
of S{u,v) (i.e e'^^"'") = S{u,v)). To treat the logarithmic discontinuities appropriately we 
define it more precisely. Let 

Hu,0 = ^^^, m>h Im^^O, (A.7) 

x{u) ^ 

and assume that u lies close to the real axis. Then the definition of o"(ti, v) is as follows: 

a{u,v) = l l-(^)+ln0(n,x(.)) Im.<0, ^^^^^ 
[ In (— a;(ti)) — ln(/>(u, Imv > 0. 

We list its most important properties below: 

a{u + ie,v) = —a{u — ie,v), |n| > 2, (A. 9) 

a{uo + i e,v) ^ lnx{v) + 2Tri@{lmv) @{uo), lio — )• ±oo, Imv ^ 0, (A. 10) 
a{-2,v)=0, a{2 + ie,v) =i7rsgn(lmv). (A. 11) 

with @{u) being the unitstep function and e is a positive infinitesimal parameter. 



Using the terminology of ref. |M| they are called canonical, simplified, hybrid etc. 
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Using the kernels K{u,v) and Kq (u — v) it is possible to define a series of kernels which 
are connected to the fermionic -functions. They are: 

KQy{u, v) = K{u - -Q, v) - K{u + -Q, v) , (A.12) 

g 9 

K^\u, v) = ^ {Kq{u -v)± KQy{u, v)) (A.13) 

and 

KyQ{u, v) = K{u, V + -Q)- K{u, V - -Q), (A.14) 

g g 

Kf{u,v) = ^(^KyQiu,v)TKQiu-v)). (A.15) 
The equation for the discontinuity function contains the kernel 

k^iu,v) = K^"'\u,v)+K^-"'\u,v), m = l,2,... (A.16) 
The kernels entering the infinite sums in the canonical equations are 

Zm au 

x{u - iO-) - x{v + i^) x{u - i^) - x{v - if) x{u + i^) 



gQM / ^ ^ ^ '9' ^" ■ ~ 9 ' ^ '9' ^" ~ 9 ' ^ ' ' g 

" x{u + i^)-x{v + i^) x{u + i^)-x{v-i^) x{u-i^) 

M-iu_^_i(Q_M + 2j) 



X 



TT (A.17) 

iJ.u-v + l{Q-M + 2j) 



and 



x(u - i^) - xiv + if) xiu + i^) - x(v + if) x(v - if) 



^-^u-v-UM -Q + 2j) 



x{u - z|) - x{v - if) x{u + i|) - x{v - if) xiv + if) 



xTT C -. (A.18) 

The equations for the momentum carrying nodes contain the dressing phase, an im- 
portant building block of the s[(2) S-matrix of the model ||50| ]. It is of the form 

^J(2)(«' ^) = S'^^'i^ - ^Qm{u, v)-^ , (A.19) 

where 2*5*^ is the improved dressing factor |||]. The corresponding s[(2) and dressing 
kernels are defined in the usual way 

KS^\i^,v) = ^^^^ogSf,^^{u,v), Kl^{u,v) = ^|;logSQM(^,^;). (A.20) 
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The source terms in the equations for the mirror magnons involve the s[(2) S-matrix 
analyticahy continued to the physical region in the first argument. 



ai(2)V ' J Sui{u-v)^,,M{n,vy 

Explicit expressions for the improved dressing factors T,qm{u,v) and T,i^m{u,v) can be 
found in section 6 of ref. [^. Their expressions contain two important functions $(xi,X2) 
and \I'(xi,X2) defined by 

^{xi,X2)=f— In— — — -, A.22 

J 2mw-X2 ^g{xi + ^ - w - ^)\ 

where the integrals run over the unit circles in anti-clockwise direction. 

The simplified equations for the momentum carrying nodes involve the kernel 

K^ = ^.-^logtQ = -KQyiIo + lQ (A.23) 

where 

oo „2 
lQ{u,v) = -Y,k2n+Q{u,v + ie) = K^^^^\u-v)-2 / dtK^^^^\u-t)K{t,v + ie), 

n=l ~^ 

K^^\u) = ^^los^-4^. (A.24) 
B. The dressing phase discontinuities 

In this appendix we present the verification of the formula ( |6.9D , which was used in the 
calculation of A in section 6. The calculation is entirely based on ref. |^], where the 
analytic continuation of the dressing phase to the mirror region and the corresponding 
integral representations were found. 

We define the dressing part of the massive node Y^^^ in the ABA limit as 



(Inyf ))drc.. = -2i e{u, Uj) . (B.l) 



'dress '^^ 

3 

This appears in ( p. 271 ). Uj are the physical particle rapidities and the second variable of 
the dressing phase 9{u,Uj) lives on the string sheet which means that the corresponding 
x^{uj) functions are evaluated in the physical (string) kinematics: 

6{u, w) = r{u, xtH) - r{u, X- H) , (B.2) 

where f"^{u,£^) is the restriction to the mirror plane of 

f{z,0 = x{x+{z),0-x{x-{z),0, ^ = xt{w), (B.3) 
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defined on the entire z torus (see Figure 1 of ^j). 

Actually, we will need this function only in three regions of the rapidity torus: 7^fc,0; 
k = 0, 1,2 We will denote the overlap of these regions with the mirror u plane with 
i?™. A; = 0, 1, 2. The regions are characterized by 

R^: Imn<-i; > 1, 



-i < Imn < |; |a;~(n)| > 1 , \x+{u)\ < 1, 
I < Imu; \x^{u)\ < 1. 



'-2 

The function xi^^O (and the complete dressing phase 9) can be expressed in terms 
of the functions ^{x,^) and ^'(x,^) defined in ref. These definitions can be found in 
appendix A. Since in our analysis the second argument of these functions (^ = xf{uj) > 1) 
plays no role in the analytic continuation process, in the rest of this appendix we will 
suppress the dependence of on ^. The function $(a;) is given by the double integral 
formula ( |A.21 ) for all |x| ^ 1 whereas for ^ix) the single integral representation ( |A.22 ) is 



valid for all |x| 7^ 1 except for an infinite number of cuts (see below). Both functions can be 
analytically extended starting from a certain region but the analytic extensions in general 
will differ from the integral formula. In this appendix $(a;),^'(x) are always understood 
as given by the integral formulae. If we analytically continue, for example, <^(x) from the 
region \x\ > 1 to values < 1, we have (for close to 1): 

^cont(x) = $(x) -^'(x), |x| (slightly) <1. (B.4) 

In the three regions we need the f{z) function, analytically continued from the physical 
7^0,0 region, is given as Q 

^>(x~), 



'^0,0 : 




= $(x+) 


7^1,0 : 




= $(x+) 


7^2,0 : 




= $(x+) 



l'(x+) + i In ^rr3|- - $(x~) + ^'(x~) . 

To calculate the discontinuity function, we also need the analytic continuation of f{z) 
from i?™ to the region i?™ through the upper cut (Imn = |) in the mirror plane. Crossing 
this cut from below, we remain in T^i^o and correspondingly we have 

Rf ■■ r{u) = ((/+).)- (n) = $(x+) - ^(x+) - cl>(x-,j , 

and since through this cut x^^^ = l/x~ this can be written as 

Ef : f*{u) = ^l^x+)-^{x+)-^{^). 

Defining (in the region just above the cut) the discontinuity 

D{u) = f{u + -)-nu + -) (B.5) 
9 9 

1 i 1 
i:>(n) = -ln-^ J -$(x) + ^'(x) + $(-), lmn>0. (B.6) 

i X — 4 X 

To obtain the discontinuity just below the real line we have to analytically continue it 
through the "slot" [—2,2]. This means that we have to continue from |x| < 1 (above the 



we find 
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real line) to |x| > 1 (below the real line). Using ( [B.4D we have 

-$(x) + ^{x) -$(x) 

and 

X X X X 

since ^'(1/a;) = ^'(a;). Just below the real line we thus have 

1 i -f 1 
D{u) = -In ^ - ^{x) +^{-) -"^{x) , lmu<0. 



(B.7) 
(B.8) 

(B.9) 



Comparing ( |B.6| ) and ( |B.S| ) we see that (since through the cut x o ^) the jump of D{u) 
through the cut is equivalent to a sign change, as expected. 

We now want to extend the discontinuity D{u) to the first Riemann sheet, i.e. the 
whole mirror plane with its infinitely many cuts at Re u > 2, Im u = Z/g, for all integers 
Z. From ( |B.6| ) we see that only ^'(m) has to be extended to the upper part of the mirror 
plane, the rest is already unambiguously defined. Similarly from ( p.9D we see that we 
have to extend —^{u) to the lower half of the mirror plane. The point is that although 
^ is already defined through its integral representation but this representation has cuts at 
the "wrong" place (|Re u\ < 2, Im u = 2Z/g, Z 7^ 0). Therefore we have to modify the 



analytic extension starting from near the real line, where D{u) is given by (B^) and (|B.S| ). 
We now perform a partial integration in ( A.22| ) and get 



iP{u) = ^'(x(u)) = - (h — ln{w -0(1-^ 



2tt 



^Onr 



dv 



InT 



1 '^Sr 



(B.IO) 



where v = w + 1/w. This can be rewritten as 

^(u) = ^{x{u)) = j> ^ H^siv) - 



■J 10 

av 



l + —(u-v) 
2 ^ ^ 



InT 



— (u — v) 
2 ^ ' 



(B.ll) 



Here the overall sign has changed since the curve 70 is defined as integration (just above 
the real line) from —2 to 2, and then back from 2 to —2 just below the real line, and this 
is a clockwise curve. From this representation we see that the cuts come from where 



2iN 



V = u- 



(first term) or 



V 



9 9 
On the upper half of the u plane the cuts are at 

2iN 



2iN , 
u-\ (second term) 



= 1,2, 



u 



Uo + 



9 



-2 < tio < 2 



(B.12) 



(B.13) 



-52- 



and come from the first term. Using the residue theorem we can calculate the jumps and 
we find that for all A'^ > 0: 



(B.14) 



2iN. 2iN ., 2iN ., 

J(uo H ) = '(p{uo H h le) — tpiuo H le) = 

9 9 9 

xiuQj — ^ 

— iln(xs{uQ + ie) — ^) + iln{xs{uQ — ie) — ^) = iln — j . 

Introducing the notation 

1 _ e 

h{u) = In (B.15) 
x[u) — 4 

we can write this jump as 

j(u) = -ih{u-^^) , lmu = —, |Re'u|<2. (B.16) 

Similarly on the lower half plane the jumps come from the second term and we have (A'^ > 0) 

2iiV, 2iN ., 2iN ., .,, , .,, 2iN , 

J[uq ) = i^luQ h le) — tpiuQ i€) = % n[uo) = i n[u H j . (B.17) 

9 9 9 9 

We can now define the modified analytic extension ip^^ (u) which has cuts at the right 
place (|Re u\ > 2) on the upper half of the mirror plane and which is defined as 

2 ( ^ 
<Imu < - : Vm (■") = ) 

2 4 r ^ 2i 

-<liau<-: il^^J^Hu) = ip(u) + ih(u ), (B.18) 

9 9 9 

4 6 2? 4i 

- < Im < - : (■") = + i K"^ ) + ^ K'^ ) ) 

and similarly for larger values of Im u. For this function we have for all iV > 0: 

[V'^^ (^)]2iV = + ^ + ^e) - V^^) (n + ^ - ie) = ih{u + ie) . (B.19) 

9 9 

Similarly on the lower part of the u plane we define ip^ (u) by 
2 

< Im u < : Vm^ (u) = ip(u) , 

--<Imu<--: ip^^u) = ip(u) + ih(u + —), (B.20) 
9 9 9 

-- <Imu < -- : i^^iu) = 'il){u) + ih{u + — ) +ih{u+ —) , 
and so on. Again, the cuts are at |Re u\ > 2 and we have for all A'' > 

u)]-2N = i'^Hu - ^ + ie) - -^^^{u - ^ - ie) = -i h{u - ie) . (B.21) 
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Recall that because of the sign change in (B^) it is —Tp{u) that has to be extended to 
the lower half plane and we have for all 7^ 



[D{u)]±2N = i h{u ± ie) = iti h{u + ie) , 



(B.22) 



i.e. all upper/lower jumps are the same and this property is also inherited by the dressing 
part of the full discontinuity function: 



[Arf(n)]±2Ar = 2j(n ± ie) 



(B.23) 



where 



. V^, xlu) x{u)-Xf^ ( Bm Rp 

J [u) = y m — i 1 — = In ' 



^(^) x{u) ^k 



R"m Bp 



In 



Y 



(0)^ 



Y 



(0) 



(B.24) 



Here = In using ( |a29D and we see from and (^]2^) that only the 



dressing part {6 part) contributes to ( B.23 ) for > 1. 

Finally we note that because of the h function parts added to V'm^ and ip'i^' these 
functions also have an infinite number of poles and zeros. We find that A^(u) has double 
zeroes at + ^^^g*"^^^ — 1)2,...) and double poles at — ^'^^^^^^ 
for all tifc. 



(iV 



-1,0,1, 



C. The asymptotic T-system: Bethe Ansatz solution 
C.l Asymptotic transfer matrices 

In the asymptotic limit the 1q,o functions tend to zero and the Y-system of AdS/CFT 
splits into two SU{2\2) Y-systems. Correspondingly two independent SU{2\2) T-systems 
generate the asymptotic solutions for the Y-functions. The asymptotic solution consistent 
with the asymptotic Bethe Ansatz equations and the multiparticle Liischer formulae 
|11] was given in |27]. 

In this appendix we discuss the form and the most important analyticity properties of 
the solution of the asymptotic T-system for states when there are A^ fundamental magnons 
with rapidities Uj present in the system. These solutions correspond to the eigenvalues 
of the fusion hierarchy of the transfer matrices built from the S-matrices of the centrally 
extended SU {2\2) algebra such that the magnon rapidities Uj play the role of the inhomo- 
geneities of the SU{2\2) vertex model. 

We introduce the following functions: 

Rmiu) = l[^^^^^^, BM = U , (C.l) 

j=i (xJP ,=1 {xJP 



RM) = U"-^^^, Bpiu) = ll^k^, (C.2) 
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where = Xs(uj ± These functions satisfy the relation 

R+{u)B+{u) = R-{u)B-{u) = (-l)^Q(n), (C.3) 

with Q{u) = rijLil^ - 

For states outside the sl(2) sector auxiUary Bethe roots appear in the formulae. To 
take into account their contribution as well, we need to introduce the following functions: 



w=n^vM^' 5K^)=n^^^, Qi{u)=Y[{u-n,), /= 1,2,3, 

j=i yyw^ j=i iyij)^ j=i 

(C.4) 

where yij = x{uij) and they satisfy the relation 

Ri{u)Bi{u) = {-lf^Qi{u), Z = 1,2,3. (C.5) 

The 3 family of Bethe roots {iti,j}i=i,2,3 correspond to the 3 levels of the SU{2\2) nested 
Bethe Ansatz. 

Using the definitions above the asymptotic solution of the T-system is given by the 
formulae as follows. 

T(°i(n) = 1, T^%) = 1, (C.6) 



^ ^ \ q[-'^] q[-'^i sir"' 



a-2 Q[a-l-2n] f bItT'^''^"'^ i?m~^~^''' \ 

+ 0(a - 2) ^ [a-2-2n] ^[a-2-2n] I n[a-2-2n] j-,[a-2-2n] ) ^^''''^ 

a-1 g[a-l-2n] / Q[a-2n] g[a-3-2n] g[a-2-2n] g[a+l-2n] \ ^ 

— Q{a — 1) ^ [a-2-2n] ^[a-2n] [ ^[a~2-2n] ^[a-l-2n] ^[a-2n] ^[a-l-2n] ) | ' 

n=0 Vs \Vl V2 ^;2 /J 

where 0(x) is the unitstep function such that 0(0) = 1. r^°^ are the eigenvalues of 
the transfer matrices corresponding to the anti-symmetric irreducible representations of 
SU (2 1 2) in the auxiliary space. The eigenvalues belonging to the symmetric representations 
are given by: 



^(0) _ ^ TT Rm I ^[-s-1] ^[s+1] Rm sr^ 

- ^f-.l^W ii r,[2j-s] W2 ^2 ^[,1 2^^^ 



where 



Vi j=i -ftp I, -ftp f.^Q 

- @is - 1) (qt'^ £ + Qt'^ Q^'' t ^.^) 

V k=0 Rp k=l / 

+ e(. - 2) Q^^l Qt'^ ^ J2 ^^^4 ' 



it=l 



g[2fc-s] g[2fc-s] 

^^''^ ~ ^[2fc-l-s] ^[2fe+l-s] ■ 
{^2 W2 



(C.9) 
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Finally T-functions on the interior boundaries of the fat-hook are given by: 



(C.IO) 





^s^=r°'l^^l (^^°^^^)' (n^l^) > ^>2, (c.ii) 

and respectively. 



where and g^°^~ are the analytic continuations of T^'^^ through the branch cut at - 

i 



They are given explicitly by: 

j7(0) _ Qs _|_ f Qt^ Q3 _|_ Q2 Qt \ _ Ql -j^g-j 



V Q2 Q2 Ql J Rp Rp Ql 

q{0) _ _Q3_ _|_ ( Qt^ Q3 I Q2 Ql \ _ Ql /q -j^^n 

Q+ B+ V Q2 Qt Q2 Ql ) BpB+Q^' 
The Bethe Ansatz equations follow from requiring that the residues of the (would-be) poles 
of the transfer matrices at the roots of the polynomials Qi^ / = 1, 2, 3 vanish: 

77+7 r^— 7 r - 1, J - 1, ••,^1 (C.15j 

Q++K,)Q,K,)Q3K,)^_^^ , = l,..,i^2 (C.16) 

Q2 (^2j) Ql (W2j) (^2j) 
Q2{U3,j) Rm{U3,j) 

The analyticity properties of the asymptotic T-functions can be easily read off from 
the formulae above. Now we summarize their most important properties. T^f^ has square 
root branch cuts along the lines Imn = ib|. J^^^^ and Q^^^ have square root branch cuts 



along the lines Imu = it^. From (C.8) it can be seen that ^ has several square root 
branch cuts between the lines Imn = but most of these cuts are generated by a gauge 
transformation and are cancelled from the Y-functions. Separating the gauge factor: 



s~l jj[2j- ^ 
= 1 -Kn 



it can be seen that f^^^ has discontinuities only along the lines Imn = 

For our considerations another important analyticity information is the large u behav- 
ior of the T-functions. Analysing the formulae above it turns out that imposing the level 
matching condition (i.e. JlfcLi i^t /^k) ~ 

rj°^(n) ~ 2(-1)'^AW ^ for |u| ^ 00 a = 1,2,..., (C.19) 
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when \lmu\ < |. The coefficient A^^^ is real and negative. In the s[(2) sector it is given 

by the simple formula = — ^ ^^t^^ , where n = ig J2 ~ is a real and positive 

number. In the region |Imtt| > | the large u behavior of T^^l{u) is different, there the 
decay is only as ^. For a = 1 we have 

rW(„) ^ ±^ 5(0) = ^(iV + /x) (C.20) 
M g 

and the upper (lower) sign is valid for u above (below) the physical strip. 

In the mirror TBA equations only the zeroes and poles of Ta,i{u) and Ti^siu) located 
in the physical strip (i.e the strip Im« < i) furthermore the zeroes and poles of T'^ and 
Q~ in the upper and lower half planes respectively are relevant. Though the asymptotic 
formulae presented above are valid for a general state in most of this paper we considered 
states where Ti^i, J^"*" and G~ have no zeros and poles in the regions |Imu| < 1/g, Imn > 
and Im-u < respectively. Wc think (but have not proved) that this condition is satisfied 
in the sl(2) sector of the theory. This set of states is definitely non-empty (see below) but 
since we are not certain if it is the entire 5l(2) sector or only a proper subset of it we have 
called it in the paper the sl{2) (sub-)sector. Here we list the properties of the transfer 
matrices that characterize this (sub-)sector and introduce some notations and definitions. 
In this (sub-)scctor there is a symmetry between the left and right wing variables so the 
wing index of the variables can be suppressed. 

In the asymptotic limit the sl{2) (sub-)sector is characterized by the properties: 

• T^^i has neither zeros nor poles in the physical strip. 

• r^°) and f^^g have no poles in the physical strip except T2^i which has poles at 
positions uj j = 1, ..,N. 

• T^'^i has AAq zeroes in the physical strip at positions ^a,j for a > 2, and we define 

i=i 

• T^^g has Ms zeroes in the physical strip at positions ^s,j for s > 2, and their contri- 
butions are taken into account by the factor fs{u) = Yl t{'^ — Cs,j)- 

• The transfer matrices have no zeroes (or poles) at the boundaries of the physical strip 
(with imaginary part zizl/g). 

• ^('^)+ has neither zeros nor poles in the region Imu > 0. 

• Q^^^~ has neither zeros nor poles in the region Imu < 0. 

N 

We also define ti{u) = Yl t{u — uj), fi{u) = l.the poles located 
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In the s[(2) (sub-)sector the source terms appearing in the TBA equations ( 3.25 - 3.31 ) 
can be expressed in terms of the Tm and fm functions as follows: 

^m|t)ui = '''mTm+2j m = 1,2,... (C.21) 

tm\w = T'mTm+2, m = 1,2, ... (C.22) 

tQ = ri^, Q = 2,3,... (C.23) 
h = 1 (C.24) 

t- = T2/f2, (C.25) 

where because of the left-right symmetry we have omitted the wing index . 

The class of minimal energy twist-two states forms an important subset within the 
s[(2) (sub-)sector. Because of the importance of these states we studied numerically the 
qualitative analyticity properties of the T-functions of this class. These can be summarized 
as follows: 

• T^^^ has neither zeros nor poles in the physical strip, 

• T2^i has 2(A^ — 2) zeros on the real axis and poles at positions uj , 

• T^^^ has 2(A^ — 2) zeros on the real axis for a > 3, 

• T^'^^ has neither zeros nor poles in the physical strip, 

• J^^^^'^ has neither zeros nor poles in the region Imu > 0, 

• Q^^^~ has neither zeros nor poles in the region Imu < 0. 

Near the small coupling limit all the magnon rapidities and the real zeroes of Ta^i are larger 
than 2, because they are of order i. 

The analyticity properties of the asymptotic T-functions summarized above determine 
the analyticity properties of the asymptotic Y- functions through the formula ( |2.3| ). Their 
state independent analyticity properties can be summarized as follows. Y^'^^ and Y^^J tend 
to constants at infinity and they have square root branch cuts along the lines Imu = 
and Imu = ib^^^i^ respectively. Y^^^ and Fg^g^ tend to —1 at infinity and have square root 
branch cuts along the lines Imu = 0, ±|, such that 1 -|- Y'/'^^ ~ l/u and 1 H — ^ ~ 1/u at 
infinity. 

The state dependent analyticity properties like the distribution of zeroes and poles can 
be read off from the local singularity structure of the corresponding T-functions studied 
above. 

C.2 Asymptotic solution for the momentum carrying nodes 

In the asymptotic limit the Y-system decouples to two SU{2\2) T-systems. Correspondingly 
we will put an L (left) or R (right) index on the asymptotic transfer matrices. Using this 
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notation the asymptotic solution for the momentum carrying nodes takes the form [27|: 

V(0) _ ( ^ 1 -T) ^3,L ^1,R ^3,R ^(0)L ^(0)R „ _ i 9 9fi^ 

N / 1,-L 3,Z^ 1,/? 3,i? 

where Je/ / = J + ^^'^ ^i.-l+^^s.-r i^i.r factor Va contains the dressing phase of 
fundamental magnons as follows. Let: 

-j N N 

B^{u) Rm{u 



with 6{u,Uj) given by (|B.2| ). Then 



a~l N 

,[a-l-2A:]^ _ 

" ■ ^ " " ■51(2)'^"'' "J^ 

fc=o i=i 



P.(n) = n Vt'-'^hu) ^llsf;iu,u,). (C.28) 



In the 5/(2) sector the following representation for Y^'^^ proved to be useful: 

y^{0) _ ^(0) rp{0)L ^io)R ^(0) _ 

'^Ifi ~ ^1,-1 -'I,! ' '^d — ydi 

where yd contains the dressing part of Y^'^q : 

+ \ J 



The Beisert-Staudacher equations follow from the requirement that ^/o|*)(^^fc) = — 1, where 

^1 ol*) ~ ((^/? )*) • words Y^^q{u) analytically continued through the branch cut at 

— I to the physical sheet and taken at the positions of the physical rapidities is equal to 

— 1. Starting from ( |C.26| ) this analytical continuation can be implemented simply by using 
( C.12| ) and ( p.l4 ). The equations take the form: 



xtV'" ^ Bi,L R3,L B^,R R3,R 
'~= ] ^^K'i) '^+~ '^+~ '^+~ "^4~ 
•^^ / -"■1,L -°3.L ^1,R ^3.R 



= -1, k = l,..,N, (C.31) 



where we introduced the notation 

AT / N 



5.i(2)(n) = n4(2)(^'^^) = n4;2)(^'^i) (c.32) 



(*) 



for short. For states outside the s[(2) sector there are 6 other Bethe equations that can be 
obtained by putting L and R indices on the 3 auxiliary Bethe equations ( C.15D -( |C.17 ). 

We close this appendix by the description the analyticity properties of Y^q. They 
have square root branch cuts along the lines Imu = ib|. They do not have any zeros in 
the complex plane and the distribution of their poles is as follows. 
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Y^'^fj has poles at positions uf^"^^ , 

Y^^Q has poles at positions vfj^^^ and vf-^^""^^^ 



For large u the Y^^q behaves as: 



a = 2,3.... 



1 



y^2( J+2+N+fj.) ■ 



for 



Imu < 



9 



The qualitative analyticity information which can be obtained from the formulae given 
in this appendix was used in the main text in the derivation of the TBA equations for the 
exact energies of these states. 

D. On the refined interpretation of the discontinuity relations 



This appendix is devoted to the detailed explanation of the refined interpretation of the 
discontinuity relations (^- 2.11 ) when local singularities lying exactly on the cut lines are 
present. We will explain the formula ( ^.12 ) in detail. To this end we will invoke the 
asymptotic solution and in order to get rid of the logarithmic type discontinuities we will 



examine the discontinuity relations for the derivative of A. We will show that (2.12) is 
valid for all states of the model and not only for states from the s[(2) (sub-)sector discussed 
mostly in the main text. 

From ( p. 26 ), which gives the asymptotic form of Y^^^ for a general state it follows that 
asymptotically A takes the form: 



aW = Ag) + a(°) + a, + A4 



(D.l) 



where 



and 



A 



(0) 
R 



A4 



(0) 



A 



"^Jeff In 2; + In 



InT, 



1,-1 



iJ+i 



B 



1,L Rz,L Bi R i?3 R 



Rl,L Rl,R B3 R 



(D.2) 



(D.3) 



This last contribution has cuts along the real axis only and thus plays no role in the 
discontinuity relations at i^y^ for 7^ 0. 



First we recall that the main result 0.23| - [B.24| ) of appendix B is that 



[A, 



d\±2N 



2 In- 



(0) 



Vln 



y 



(a)(0) 



y 



(a)(0) • 



(D.4) 



Using ( p. 71) of appendix C we will show that if we forget about the local singularities 
of the asymptotic T- and y-functions we get the naive equality for the derivatives: 



A 



(0)/ 
R 



±2N 



D 



{±2N){0)i 
12R ' 



(D.5) 
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where 



D 



(±2JV)(0) 
12B. 





( 




± 













(+)(0) 



Y- 



+ 



±2N 



N 

E 

771=1 



In 1 + 



(+)(o) 



Y 



(D.6) 



and similarly we get analogous expressions for This complicated expression of Y- 

functions becomes much simpler in the language of T- functions. For the sake of simplicity 
let us consider the discontinuity relation in the upper half plane. Similar considerations 
that lead to (6.11) in the main text give 
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Using ( C.12| ) we write 
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From the expressions of appendix C it can be seen that rj^m^ j^g^g 

no square- 
root discontinuities in the upper half plane and in the strip |Imii| < y respectively, and it 
follows that 
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which implies ( D.5| ). This together with ( p.4[ ) would justify formula ( |2.9D for the disconti- 
nuities of A. 

Now let us examine whether local singularities of T-functions lying exactly on the 
cut lines might modify ( |2.9D ? To investigate this point some information on the local 
singularities of the T-functions is necessary. Since we require the y-functions to be real 
analytic functions of u their local singularities can be either real or come in complex 
conjugate pairs. Due to the T-representation of the y-functions, the local singularities of 
the T-functions generate this structure. Thus, T-functions have the same local singularity 
structure apart from some irrelevant modification coming from the fact that the asymptotic 
form of the T-functions given in appendix C corresponds to a gauge where the T-functions 
are not real analytic. 

From this discussion it is clear that the real zeroes are the most important objects for 
the question we are investigating here since they preserve their imaginary parts, while the 
other singularities move in the complex plane if we tune the coupling g from low to higher 
values. In principle it can happen that at certain values of g some complex roots lie exactly 
on one of the cut lines, but since their imaginary parts are not restricted by any symmetry 
they move off the cut line if g is changed a little. 

For the sake of simplicity we exclude this zero measure set of the possible values of g 
from our discussion and then we have to consider whether the real zeroes of the T-functions 
modify (|2.9| ). Formula ( p. 7)) implies that the real zeroes of the T-functions can lie on the 
cut lines only in the term In' T^^^] -, and only when = 1 . This means that 
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where the set {?^2^''}^i denotes the real zeroes of T^^^' with absolute value larger than 2. 
The equality can be restored if we remove the contribution of these local singularities from 



(+) 



RiO) 



the right hand side. For this purpose we introduce the polynomial (^) = " 

k=l 

and the modified equality which holds for all |n| > 2 takes the form: 
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In order to rephrase this modification in the language of gauge invariant 1"- functions, let 
us recall the naive discontinuity relation ( p.SD for = 1: 
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Since 1 + 



y(+) 
1 \ vw 



^21^21' clear that this term is responsible for the neccessary 



modifications. On the other hand Y^^' ~ T2^i, thus the zeroes of ^2^^ can also be defined 
as the zeroes of y^"*"^ lying on the real cut line. 

Similar considerations can be applied for the discontinuity relations corresponding 
to the lower half plane and to the left part A^. Putting together the results of these 
considerations it turns out that the form of the proposal ( |2.9| ) is unaffected by the local 
singularities for N >2, but must be modified when A'^ = 1 as follows: 
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where P2'\u) is the polynomial having zeroes at the positions of the real zeroes of y^"^ 
with absolute values larger than 2. 

Similar arguments as we applied for A shows that the other discontinuity relations 
( |2.10 ) and ( 2.11| ) are not to be modified. We note that the problem of local singularities 
lying on the cut lines is not only of academic interest, but also occurs in realistic cases, for 
example for the minimal energy twist-two states in the s[(2) sector of the model. 

A final remark on the modification of ( ^.9| ) is in order. If in (|2.9| ) the term discontinuity 
is not defined by , but is defined as that part of the discontinuity which comes from the 
square-root behavior of the functions, then (|2.9D is correct for = 1, too. The subtlety is 
associated with the log type branch points lying on the cut lines which becomes important 



when (2.9) is translated to dispersion relations. 



E. Regularity and the exact Bethe equation 

In this appendix we clarify the consequences of imposing the exact Bethe equations and 
compare this to what is known in analogous relativistic integrable models. 

In simple relativistic integrable models we have only one massive node (with corre- 
sponding Y-system function denoted Yq) and, like for the Sine-Gordon model, the XXX 
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model, or (with some small modification) for the 0(4) nonlinear sigma model, the only 
Y-system equation containing the factor 1 + Yq is for the neighboring Y-function Yi and is 
of the form 

Y+Y^- = {l + Yo){l + Y2). (E.l) 

For the asymptotic Y-system we have to drop the (l+lo) factor since Yq = asymptotically. 

The zeroes of Yi are the physical rapidities {Oj} and the elements of the set {r^}. Both 
sides of ( [E.lD vanish at the points r^. For the right hand side of the equation this vanishing 
is due to the factor (l+l2)- The above conclusion is the same whether we consider the exact 
problem or the asymptotic limit. In contrast, the behavior of the two sides of the equation 
at singular points associated to the physical rapidities 6j is very different asymptotically 
and for the exact solution. Asymptotically the right hand side is regular at 9^ and since 

[±2l 

one of the factors vanishes on the left hand side, the other factor has to have a pole at 9j 
to compensate. However, for the exact solution the analog of ( [).12D 

1 + Yoief) = (E.2) 

is imposed and now both sides vanish at 9^. There is no need for compensation and Yi 

remains regular at ^^^'^^ • Thus as a consequence of imposing the exact Bethe equation ( [E.2[ ) 
the neighboring Y-function becomes more regular. This occurs two steps away from the 
real line. 

The consequences of the exact Bethe equations (9.12) are very similar, but we have 



to take into account that in the AdS/CFT case these quantization conditions are imposed 
after analytical continuation to a different sheet, the string sheet. 



The relevant equation is ( 3.25 ). We rewrite its inverse after the analytic continuation: 



^ = (l + yi(,))A^(,), (E.3) 



where 



1 + y/"^ 

M = (E.4) 



~ l\vw 

and we recall that for any function / we have 

/m = ((/")*)+• (E.5) 



Again, as for the relativistic case, the first factor on the right hand side of ( |E.3| ) is absent 
in the asymptotic limit. 

Using the explicit formulae given in appendix C, we can convince ourselves that the 
factor A^{*) is regular at singular points associated to the physical rapidities Uj. We also see 
that the first factor on the left hand side of (|E.3| ) has zeroes at {uj}. (For the asymptotic 
case this can be seen using the explicit formulae of appendix C again and for the exact 
problem it follows from the discussion of subsection 5.2). 

The consequences of ( |E.3D are very similar to what we have seen above. Asymptotically 
^(^) has to have poles at uj to compensate the zeroes of the first factor. In the exact case 
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because of ( 9.12 ) the right hand side also vanishes and ^(*) remains regular. Thus also 
in our case the exact solution is more regular than the asymptotic limit. We note that 
this occurs at two TBA steps away from the real line in the physical (string) sheet: while 
asymptotically (y1"^ )* has zeroes at uf, these zeroes are absent for the exact solution. 
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